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(lhaptof  1 
INIKODIICTiON 


1.1  liil  ion 

Tlio  I'urposo  o)  Mils  report  is  to  jnovide  an  cn{j i neer i nj;  basis 
for  pred  ict  iiit;  I  lie  elfccts  of  nearby  power  transmission  lines  on  the 
l)orfoi mam  ('  of  Noii-li i  reef  i ona  I  Itcacuns.  Iwo  l)asic  mechanisms  for 
possible  I  lit  erforeni  (’  liave  been  considered;  Ibulio  l're»iuency  noise 
goiier.i  i  ed  by  t  lie  lines  lliemselvos,  ant!  passive  reradialion  of  the  NDH 
si(;iial  by  I  he  t  r.insiii  i  -s  ion  line  lowers  and  conductors.  Since  the 
proiluciion  of  radio  noise  is  very  different  for  Ai;  and  l)b  lines,  a 
separaie  chapter  is  devoted  to  each.  A  third  chapter  is  devoted  to 
lerad  i.il  u'li  effects  (both  AC  and  DC  linos  will  behave  similarly). 

■file  cone  111.  ions  re.ichocl  tire  discussed  in  a  final  chapter;  however, 
the  I  undainenta  1  eoiulusion  is  that  locating  the  NDU  triinsmitter  near 
a  poser  transmission  line  should  not  have  atiy  detrimental  effects  due 
to  tin  alien  (•  mecha  ii  i  sms .  It  must  be  immediately  pointed  out,  however, 
that  Ollier  mcch.ini sms ,  such  as  jiowerline  carrier  radiation,  arc  not 
included  in  the  above  statement. 

In  this  report  the  siyiial  from  the  NUB  transmitter  is  the  de¬ 
sired  signal,  and  the  Kl  noise  from  the  powerlines  or  the  reradiuted 
signal  from  the  powcrline  structures  the  undesired.  The  FAA  Handbook^ 
states  that  if  the  level  of  the  undesired  signal  is  15  dB  below  that 
of  the  desired  signal,  the  AUf  positioning  error  will  be  less  than  ±l 
degree.  This  criteria  will  be  applied  throughout  this  work  for 
assessing  the  potential  of  interference. 


The  objective  of  this  report  is  to  provide  a  theoretical  capabi¬ 
lity  to  predict  the  critical  distance  between  the  powerlines  and 
receiving  aircraft  where  the  ratio  of  the  desired  signal /undesircd 
noise  is  15  dB  as  a  function  of  effective  radiated  power  of  the  NOB 
transmitter,  relative  permittivity  and  ground  conductivity  of  the  earth, 
powerline  parameters,  distance  separating  the  NDB  transmitter  and  the 
powerline,  evaluation  of  the  ADF  receiver  and  the  NDB  frequency.  This 
objective  has  been  successfully  met.  In  particular,  mathematical 
models  now  exist  which  will  provide  the  desired  predictions. 


Chapter  II 

PREDICTION  OF  RADIO  INTERFERENCE  NOISE  FROM  AC  POWERLINES 


3 


2.1  Introduction 

Radio  Interference  (RI)  noise  from  AC  powerlines  originates 
from  corona,  which  is  defined  as  "a  luminous  discharge  due  to  ioni¬ 
zation  of  the  air  surrounding  a  conductor  around  which  exists  a 

2 

voltage  gradient  exceeding  a  certain  critical  value".  An  IEEE 
Report^  states  that  the  Rl  performance  can  be  divided  into  three 
stages;  1)  generation,  2]  propagation  and  3)  radiation. 

The  generation  of  corona,  and  therefore  RI,  is  basically  a 
function  of  factors  which  can  be  divided  into  two  categories: 
i)  lino  design  factors  and  ii)  atmospheric  and  environmental 
factors.  The  former  includes  the  diameter  and  spacing  of  conduc¬ 
tors,  phase  spacing  and  phase  configuration.  The  latter  deals 
witli  the  accumulation  of  foreign  particles  on  the  line  conductors 
and  weather  conditions. 

'Ihe  propagation  of  RI  along  the  powerlines  is  primarily  a 
function  of  the  line  design  and  earth  resistivity,  both  of  which 
influence  the  rate  of  attenuation  of  the  traveling  waves.  The 
radiation  of  RI  away  from  the  powerlines  is  essentially  dependent 
on  the  positions  of  th’  phases  of  the  lines  with  respect  to  the 

location  of  the  observer. 

2  3  4 

Studies  and  reports  ’  ’  have  consistently  concluded  that  the  RI 


noise  level  is  at  its  peak  under  heavy  rain  condition.  All  predictions 
of  RI  noise  from  AC  powerlines  given  in  this  report  will  be  for  the 
worst  case  condition  of  a  heavy  rain  shower. 
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2.2  Computation  of  RI  Noise  from  AC  Powerlines 

4 

CIGRn/lEEE  SuT*vey  Results  list  various  methods  for  the  calcu¬ 
lation  of  RI  noise  from  AC  powerlines.  They  are  grouped  into  two 
categories:  i)  the  comparative  method  and  ii)  the  analytical  method. 

The  former  makes  use  of  an  empirical  formula  developed  after  a  careful 
measurement  study  with  particular  reference  criteria.  To  predict  the 
RI  noise  level  from  lines  of  different  designs,  various  corrections 
for  corona  generation,  measurement  frequency  and  lateral  distance  are 
made  based  on  measurements  of  the  variables  involved.  The  analytical 
method,  on  the  other  hand,  makes  use  of  a  theoretical  characteristic 
quantity  of  RI  generation  which  is  called  the  excitation  function,  and 
proceeds  to  compute  the  total  corona  currents  on  the  lines  and  the 
resulting  RI  noise  field  strength  at  the  location  of  the  observer. 

The  method  that  will  be  used  here  belongs  to  the  latter  group 

5 

and  is  called  the  "modal  analysis  method".  It  assumes  that  the  power- 
lines  are  a  reflection-free  system  of  n  ideally  parallel  conductors 
at  constant  height  above  a  perfectly  flat  conducting  ground  whose 
potential  is  zero.  The  influence  of  the  ground  wires  is  included 

in  determining  the  surface  gradients,  but  is  neglected  for  the  modal 
currents.  The  method  is  composed  essentially  of  the  following  steps: 

Step  1.  Determination  of  the  maximum  bundle  surface  gradient 
The  maximum  bundle  surface  gradient  for  each  phase  can  be  determined 
by  using  a  method  developed  by  Markt  and  Mangel e^.  The  calculation 
is  carried  out  with  the  assumption  that: 

the  ground  is  an  infinite  horizontal  conducting  plane  surface, 
the  conductors  are  smooth  infinitely  long  circular  cylinders 
paiillcl  to  each  other  and  to  the  ground  plane. 
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the  conductors  are  equipotential  surfaces,  with  known  potentials 
applied  to  them,  the  ground  plane  is  assumed  to  be  at  zero 
potential , 

the  influence  of  powerline  towers  and  of  any  objects  in  the 
vicinity  is  neglected  and 

the  horizontal  spacing  between  the  conductors  remains  constant 
and  the  height  of  the  conductors  above  ground  is  also  constant. 
The  procedures  for  the  calculation  of  the  maximum  bundle  surface 
gradient  are  detailed  below: 

Step  a:  Hach  conductor  bundle  is  replaced  by  an  equivalent  conductor 


having  a  radius  r  defined  by 
eq 


r  =  (n  a  a""^) 
eq 


1/n 


(2-1) 


where :n  is  the  number  of  subconductors, 
a  is  the  subconductor  radius, 

A  is  the  bundle  radius. 

Step  h:  With  the  bundles  represented  by  the  equivalent  conductors, 

the  total  charge  on  each  of  them  is  calculated  by  the  Maxwell 
potential  coefficient  method  assuming  appropriate  potentials 
on  the  different  phases.  Any  ground  wires  are  also  taken  into 
account.  Program  CHARGE  listed  in  Appendix  B  is  used  to 
compute  the  values  of  the  total  charge  on  the  bundles. 

Step  c:  The  average  bundle  gradient  is  calculated  as: 


E 


‘It 


av  2TrG 


1_ 

na 


(2-2) 


Step  d:  I'inally  the  maximun  bundle  surface  gradient  is  obtained  as: 

.  (-1)  t)  (.’-3) 


0 


■Step  2.  Ucterminat  ion  of  lU  excitation  function  under  heavy  rain 
conditions 

The  excitation  function  is  considered  to  be  the  specific  measure  of 
the  cause  of  KI.  It  is  related  to  the  induced  corona  currents  by  the 
equation^ 


where:  i  is  the  induced  corona  current  injected  per  unit  length  of  the 
conductor, 

C  is  the  capacitance  of  the  conductor 
is  the  permittivity  of  air, 
r  is  the  excitation  function. 

g 

In  this  analysis,  an  empirical  formula  developed  by  EPRI  will  be 
used.  For  each  phase,  the  value  of  excitation  function  1'  is  given  by: 

r(n,d)  =  78.0  -  S80.0  38.0  log  (  d  )  +  K  (2-5) 

E  3.8  " 

m 

where  r(n,d)  is  the  excitation  function  for  a  bundle  of  n  subconductors 
of  diameter  a  (cm)  in  units  of  dB  above  1  UA/m**, 


K  =  7  dB,  if  n  =  1, 
n 

K  =  2  dB,  if  n  =  2, 
n 

K  =  0  dB,  if  n  ^  3. 
n 

Step  3.  Uetermination  of  the  geometric  matrix  [G]  of  the  powerlines 
From  the  powerline  configuration,  the  elements  of  the  geometric  matrix 

g 

[Gj  can  be  calculated  with  the  equations 


2h. 

1 

r . 

(2-6) 

1 

I).  . 

_JLL 

d 

(2-7) 

ij 
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where:  h.  is  the  height  above  ground  of  the  ith  conductor, 

rj^  is  the  radius  of  the  ith  bundle  conductor, 

is  the  distcince  between  the  ith  conductor  and  the  image  of 

the  jth  conductor  in  the  ground  plane, 

d. .  is  the  distance  between  the  ith  and  jth  conductors, 
ij 

Step  4.  Determination  of  the  modal  transformation  matrix  [M] 

The  U1  propagation  along  the  powerlincs,  including  the  corona  generated 
current  along  the  conductors,  can  be  written  by  the  equations^ 


[VI  =  -  [z][Il 


(2-8) 


[Tl  =  -  fyim  >  [J] 


(2-9) 


where:  [V]  is  the  column  vector  of  voltage  on  the  line, 

[I]  is  the  column  vector  of  current  on  the  line, 

[J]  is  the  corona  current  density  injected  on  the  conductor, 

[z]  is  the  square  matrix  of  series  impedence  per  unit  length  of 

the  line, 

[y]  is  the  square  matrix  of  the  shunt  admittapfs.  per  unit  length 
of  the  line. 

Modal  analysis  is  used  to  simplify  the  equations  (2-8)  and  (2-9) 
above  into  a  number  of  uncoupled  sets  of  equations  which  can  be  solved 
conveniently.  At  this  stage,  it  is  assumed  that  the  line  is  lossless. 
The  lossless  line  parameters  are  given  by 

[7-1  =  a)[hl  =  u)  ^  [r.]  f2-10) 

[y]  =  (Die]  =  u)  2ire^lG]‘^  (2-11) 

where:  [L]  is  the  inductance  matrix  of  the  line, 

[C]  Ls  the  capacitance  matrix  of  the  line. 
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is  the  permeability  of  air, 
is  the  permittivity  of  air. 

Equation  (2-4)  in  Step  2  can  be  rewritten  as 

f-'l  =  2^^  (2-12) 

o 

Substituting  equations  (2-10)  to  (2-12)  into  equations  (2-8)  and  (2-9) 
yield 

jfm  .-5^(G1[I1  (2-13) 

(I) .  -  u2«E_jinr‘i''i  *  (or'tn  c^-w 

Let  the  modal  matrix  of  [G]  be  [M]  and  [X]j  be  the  diagonal  spectral 
matrix  of  [C] ,  then 

fM]’^[G][M]  =  [X]^  (2-15) 

or 

(2-16) 

or 

[[(-,]-  [X]^}  [M]  =  0  (2-17) 

The  equation  (2-17)  above  has  a  non-trivial  solutions  only  if  the 
determinant  of  {i«i  -  is  zero.  Solving  this  characteristic 

equation  of  [G]  yields  the  eigenvalues  X^,  X^  and  X^.  Substituting 
these  back  into  equation  (2-17)  will  yield  corresponding  eigenvectors. 
Normalized  column  eigenvectors  will  form  the  modal  transformation 
matrix  [M] . 

Step  5.  Determination  o f  modal  components  of  corona  current  densities 
in  the  conductors 
i.et ; 

fVl  ^  [M]fVJ 


.  _ 


(2-18) 
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[1] 

(2-19) 

I'l 

=  fM][J^] 

(2-20) 

fn 

(2-21) 

where  V^,  and  are  the  modal  components  of  the  voltage  and  currents 
and  is  the  modal  components  of  the  RI  excitation  function.  Substitii- 
ing  equations  (2-18)  to  (2-21)  into  equations  (2-13)  and  (2-14)  yield 

J7  [n][Ml[T^]  (2-22) 

[MUTJ  =  -  a,2Tie^[G]'l[M][V^]  ^  [R]'^[M][rj  (2-2.3) 

Since  UT  is  a  direct  result  of  the  induced  corona  currents  in  the 
conductors,  only  the  second  part  of  equation  (2-23)  will  be  a  matter 
of  interest.  The  modal  components  of  induced  corona  currents  are 
therefore  obtained  by  the  relationship 


f.ij  =  rM]-^G]’^(M]rrj 

(2-24) 

(jj  =  [M]'^G]'^[r] 

(2-25) 

With  the  known  values  of  the  elements  of  matrices  [M]  and  [G]  and 
the  excitation  function  for  each  phase,  the  values  of  the  modal  com¬ 
ponents  of  corona  currents  for  each  phase  can  be  determined. 

Step  6.  Determination  of  corresponding  modal  components  of  the  currents 
in  the  conductors 

At  this  point,  the  effect  of  losses  will  be  introduced  in  the  form 
of  modal  attenuation  factor  a  which  is  a  function  of  frequency  and 
ground  resistivity.  Required  values  of  the  attenuation  factor  can  be 

g 

obtained  by  using  the  empirical  equation 


10 


=  (f”-®  )  a'"(l. 0,100)  (2-26) 

where :m  is  the  mode  number, 

f  is  the  frequency  in  Megahertz, 

P  is  tlic  ground  resistivity  in  Ohm  m 

8  m 

From  the  HPRI  Reference  Book  ,  the  values  of  a  (1.0,100)  are  reproduced 
in  Table  2.1  shown  below. 


Voltage 

class 

1 

1 

Attenuation  Constants 

7  -1 

1 

(kV) 

a 

(m' 

) 

a  (m 

) 

a 

(m' 

) 

362 

8.0 

X 

10-" 

60.0  X 

lo-'^ 

350.0 

X 

lo-"^ 

550 

9.3 

X 

10'^’ 

70.0  X 

10-^’ 

350.0 

X 

10-^’ 

800 

10.0 

X 

10'^’ 

70.0  X 

10-^^ 

350.0 

X 

10-" 

1200 

10.6 

X 

10’^ 

84.0  X 

10-" 

350.0 

X 

10'^ 

1500 

10.6 

X 

10-^’ 

84.0  X 

10-" 

350.0 

X 

10-^ 

Table  2.1  Modal  attenuation  constants  for  AC 
powerlines  at  1  MMz  with  ground 
resistivity  being  100  Ohm-m. 

Under  heavy  rain  conditions,  the  value  of  p  is  assumed  to  be  75.00 


Ohm-m.  Values  of  the  modal  components  of  the  currents  in  the 
conductors  arc  determined  by  the  equation*'’ 


[I 


(2-27) 


Step  7.  Determination  of  the  corresponding  phase  currents 
Refer ing  to  equation  (2-19)  in  Step  5,  it  can  be  rewritten  as 

fT]  =  [MlfT^"’]  (2-2R) 

or  it  ean  be  expanded  to  become 
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”l2 

”l3 

^1 

CM 

<M 

”23 

.^31 

^32 

”33 

n 


(2.29) 


Equation  (2-29)  indicates  that  each  modal  current  component  flows 

2 

in  all  the  three  phases.  For  example,  modal  current  will  result 
2  2 

current  M,_T  flowing  in  conductor  1,  in  conductor  2  and 

12  c  22  c 

2 

in  conductor  .1.  The  current  components  in  the  three  conductors 
of  any  given  mode  are  in  phase  with  each  other.  However,  the  currents 
corresponding  to  different  modes  are  not  in  phase  due  to  different 
velocities  of  propagation. 


Step  8.  Determination  of  the  R1  noise  field  strength  with  the  location 
of  the  observer  at  ground  level 

With  the  known  values  of  the  corona  currents  in  each  phase  computed 
above,  the  resulting  magnetic  field  and  electric  field  can  be 
conveniently  calculated.  Let  the  position  of  the  observer  be  at  P 
as  illustrated  in  Fig.  2.1.  The  magnetic  field  strength  at  this  loca¬ 
tion  with  respect  to  the  center  phase  will  be 


2h 


u2  ^  2 

h  +  X 


(2-30) 


Assuming  a  TEM  mode  of  wave  propagation,  the  corresponding  electric 


field  strength  is  obtained  by  the  relationship  of  E  =  Z^ll,  where 


is  the  wave  impedance  of  free  space.  Since  Z^  =  IZOtt  ,  the  electric 


field  strength  then  is 
...  2h 


S 
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where  f)0*  — :5 - TT  is  termed  as  the  field  factor. 

h^^  + 

Let  the  field  factors  for  conductors  1,2  and  3  be  Fj(x),  F^Cx)  and 
F^(x)  respectively.  Refering  to  Fig.  2.1  and  equation  (2-31),  the 
field  factors  for  respective  conductors  will  be 


,  120h 

^  h^  4  (x4D)'‘ 

(2-32) 

,,  ,  ,  120h 

=  ,2  2 

(2-33) 

h  4  X 

..  ,  ^  120h 

h^  4  (x-D)^ 

(2-34) 

Using  the  above  equations  (2-32)  to  (2-34),  the  electric  field  strength 
at  any  point  P  due  to  corona  for  each  mode  can  be  obtained  by 

=  Fj(x)  .1™  4  r2(x).i®  4  F3(x).t™  (2-3.-;) 

where  m  is  the  mode  number  1,2  or  3.  The  total  electric  field  strength 
for  each  phase  is  calculated  from 


Finally,  the  resultant  electric  field  strength,  therefore  RI  noise, 
due  to  corona  currents  on  all  the  phases  at  any  point  P  is 

(2-37) 

where  F  is  in  pV/m,  or  it  can  be  converted  into  dB  above  1  pV/m. 
Appendix  A  outlines  tbe  derivation  for  equations  of  field  factors  for 
point  P  anywhere  in  space  and  Appendix  C  details  the  Program  ACRI 
listings  for  the  computation  of  RI  noise  level  from  AC  powerlines. 


The  RI  noise  level  computed  here  is  for  1  kHz  receiver  bandwidth. 

In  general,  a  bandwidth  correction  factor  of  10  log  b(dB),  where  b 
is  the  receiver  bandwidth,  is  added  to  the  result  obtained  in  equation 
(Z-."??).  This  form  is  applicable  throughout  this  report.  Table  2.2 
lists  typical  values  of  bandwidth  correction  factors. 

b  (kHz)  10  log  b(dB) 

1  0.00 

2  3.01 

.•?  4.77 

4  6.02 

5  6.99 

Table  2.2  Typical  values  of  bandwidth  correction  factor 
For  illustration,  some  actual  line  designs  are  considered.  Table  2.7> 
lists  the  line  voltages  and  parameters  and  Fig.  2.2  show  the  line  con- 
f igurat ions 


Line 

Conf ig. 

Line 

Voltage 

(kV) 

H  (m) 

D  (m) 

d^  (cm) 

n 

A  (cm) 

«  1 

.345 

13.61 

8.31 

3.038 

2 

45.70 

«  2 

500 

14.43 

12.19 

2.959 

3 

52.80 

#  3 

765 

20.83 

13.72 

2.959 

4 

64.70 

n  4 

1100 

21.. 34 

15.24 

3.556 

8 

101.60 

Table  2.3  AC  line  voltages  and  parameters  considered 

where  II  is  the  average  height  of  the  conductors, 

0  is  the  spacing  between  the  conductors, 

cl  is  the  diameter  of  the  subconductors, 

c 

n  is  the  number  of  .subconductors  in  a  bundle. 


I 


«.* ... 
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A  is  the  diameter  of  a  bundle. 

By  using  the  modal  analysis  method  detailed  above,  the  RI  noise  level 
is  calculated  for  each  of  the  line  designs  listed  in  Table  2.3  at 
frequencies  of  200  kHz  and  500  kllz  with  the  observer  at  various  alti¬ 
tudes.  Fig.  2. ,3  to  2.6  illustrate  the  characteristics  of  the  calculated 
UT  noise  levels  radiated  from  the  AC  powcrlines. 

2.3  Computation  of  Signal  Strength  from  an  NDB  Transmitter 


Based  on  the  assumptions  that  the  heights  of  the  transmitting  and  re¬ 
ceiving  antenna  arc  within  a  few  wavelengths  of  the  earth's  surface 
•and  that  they  arc  separated  by  a  short  distance  such  that  the  curvature 
of  the  earth  is  negligible,  a  flat  earth  model  can  be  used  for  the 
computation  of  the  electric  field  strength  of  the  signal  from  an  NDB 
transmi ttcr . 

An  initial  approximate  computation  as.sumes  that  the  earth  is  infinitely 
conducting.  Later,  this  result  is  modified  by  multiplication  by  the 

9 

flat  earth  attenuation  function  developed  by  Wait  ,  which  will  give  tlic 
final  result  of  the  electric  field  strength  over  a  finitely  conducting 
earth.  Fig.  2.7  illustrates  a  vertical  electric  monopole  of  length  "t" 
carrying  current  I  located  on  the  surface  of  a  finitely  conducting 
earth. 

Ilarrington^^  states  that  far  from  a  current  element  in  free  space,  the 
electric  field  strength  is  given  by 


E  =  n  •  si 

6  ^  2Ar  ° 


sin  0 


(2-38) 


and  the  actual  average  power  radiated  by  the  current  element  is 

lu 


"r  ■ 


2 


(2.39) 
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Fig,  2.3a  Calculated  RI  noise  profile  for  345  kV  AC  pouerline  at  200  kHz 
under  heavy  rain  condition 
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under  heavy  rain  condition 
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under  heavy  rain  condition 


2.5b  Calculated  R1  noise  profile  for  765  kV  AC  powerline  at  500  kHz 
under  heavy  rain  condition 


Fig.  2.6a  Calculated  RI  noise  profile  for  1100  kV  AC  powc'-line  at  200  kHz 
under  heavy  rain  condition. 


where:  n  is  the  wave  impedance  of  free  space, 

X  is  the  wavelength, 
k  is  the  wave  constant 

r  is  the  distance  from  the  transmitting  antenna. 

Since  we  are  interested  only  in  the  magnitude  of  the  electric  field 
strength  near  the  ground,  equation  (2-38)  can  be  simplified  into 


If  the  earth  is  assumed  to  bo  present  but  infinitely  conducting,  due 
to  image  theory,  the  far  field  term  will  be  exactly  doubled.  Therefore 

i;„  .  n  II  (2-41) 


and  the  corresponding  value  of  the  actual  average  power  radiated  is 
given  by 


(2-42) 


From  equation  (2-42)  above. 


j4niT 


Substituting  equation  (2-43)  into  (2-41)  yields 


(2-43) 


r  =  9  487 

^'0  Jatj  r  r 


(2-44) 


Finally,  by  multiplying  equation  (2-44)  by  the  flat  earth  attenuation 
function,  the  resulting  electric  field  strength  over  a  finitely  con¬ 
ducting  earth  is 


P  =  - —  (1.00  -  6c*  erfc(z)) 


(2-4.^) 


where:  F.  is  the  electric  field  strength  in  yV/m, 
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Is  the  effective  radiated  power  in  watts. 


<S 

and  n  =  e 

r 


Jn  -  1 
n 

.  1.8  X  10^*  a 

3  - T - 


(2-46) 


(2-47) 


with  and  o  the  relative  permittivity  and  conductivity  in  mho/m  of 
the  earth  respectively,  f  the  frequency  in  kilohertz, 


and  I)  =Jr^  +  (h^  -  h^)^ 


(2-48) 

(2-49) 


with  k  =  2ti/X,  X  being  the  free  space  wavelength,  and  h^  and  h2  the 


antenna  heights;  erfe  is  the  complementary  error  function  as  defined 


11 


by  Abramowitz  and  Stegun  and 

j'rr/4 
z  =  e-^  • 


^1  "  *^2 

6  .  (i.oo.-i^) 


(2-50) 


Assuming  that  the  transmitting  antenna  is  a  40-foot  high  vertical  radia¬ 
tor,  the  relative  permittivity  of  earth  is  10.0  and  the  ground  conducti¬ 
vity  he  10.0  mmho/m.  Fig.  2.8  to  2.11  illustrate  the  patterns  of  the 
electric  field  strength  for  NDB  transmitters  with  an  effective  radiated 
power  (r:RP)  of  0.05  watt,  0.10  watt,  0.5  watts  and  1.0  watt  re.spective- 
ly,  at  the  frequencies  of  200  kHz  and  500  k!Iz  with  the  receiver  altitude 
at  ground  level  and  at  1500  feet  above  ground.  Appendix  D  details 
program  signal  listing  for  the  computation  of  the  electric  field 
strength  from  an  NDB  for  short  distances  under  the  above  conditions. 


2.4  Prediction  of  the  Critical  Distance  Where  the  Ratio  of  Desired 
Signal/Undesircd  Noise  is  15  dB 

The  desired  signal  is  the  signal  from  the  NDB  transmitter  and  the 


200  kHz 


Transmitter;  ERP  =  0.05  watt,  e  =  10.0,  ground  a  =  10.0  mmho/m 
f  =  500  kHz.  ^ 


500  kHz 


Fig.  2.10b  Transmitter;  ERP  =  0.5  watt,  e  =  10.0,  ground  a  =  10.0  mmho/m 
f  =  500  kHz. 


Fig.  2.11b  Transmitter;  ERP  =1.0  watt,  e  =  10.0,  ground  a  =  10.0  mmho/m 
f  =  500  kHz. 
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undesired  signal  is  the  RI  noise  from  the  AC  powerlines.  The  former 
has  been  dealt  with  in  Section  2.3  and  the  latter  in  Section  2.2. 

Consider  the  scenario  illustrated  in  Fig.  2.12.  The  scenario  in 
question  illustrates  a  worst  case  condition  where  the  path  of  the 
aircraft  flying  over  the  powerline  is  perpendicular  to  it,  the  air¬ 
craft  and  the  NDB  transmitter  are  on  the  opposite  side  of  the  power¬ 
line  and  a  heavy  rain  is  falling.  Let  the  distance  separating  the 
NDB  transmitting  antenna  and  the  powerlines  be  y  nautical  miles,  the 
lateral  distance  between  the  powerlines  and  the  location  of  the 
receiver  (aircraft)  be  x  feet  and  let  the  receiver  be  at  h  feet  above 
ground.  For  various  line  designs  mentioned  in  Section  2.2  and 
different  values  of  ERP  of  NDB  transmitter  covered  in  Section  2.3, 

Fig.  2.13  to  2.24  indicate  the  critical  distance  x  feet  where  the 
ratio  of  the  desired  signal /undesired  noise  is  15  dB  with  the  receiver 
at  various  altitudes. 

For  example,  in  Fig.  2.23a,  where  the  ERP  of  the  NDB  is  1  watt  and  the 
line  voltage  is  765  kV,  if  the  distance  separating  the  NDB  and  the 
powerline  is  10  NM,  the  critical  distance  for  an  aircraft  flying  at  an 
altitude  of  1500  feet  is  about  300  feet  from  the  line. 

After  examination  of  Figures  2.13  to  2.24  it  will  be  noticed  that  the 
critical  distance,  and  thus  the  RI  noise  level,  increases  as  line 
voltage  is  increased  from  345  kV  to  500  kV  to  765  kV,  but  decreases 
as  the  voltage  is  further  increased  to  1100  kV.  The  reason  for  this 
is  that  the  1100  kV  line  considered  has  8  subconductors,  while  the 
765  kV  line  has  only  4  (see  Fig.  2.2,  Table  2.3).  This  causes  the 
maximum  electric  field  for  the  1100  kV  line  to  be  less  than  that  for 
the  765  k\'  line,  thus  reducing  the  RI  noise  level. 
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2.5  Conclusion 

At  ground  level,  the  radiated  RI  noise  is  vertically  polarized.  How¬ 
ever,  at  any  other  location  in  space,  the  radiated  RI  noise  is  not 
entirely  vertically  polarized.  In  determining  the  critical  distance, 
which  involves  the  computation  of  desired  signal /undesired  noise  ratio, 
no  attempt  is  made  to  take  the  RI  noise  component  aligned  with  the 
vertical  direction.  Instead,  a  worst  case  condition  is  considered  in 
which  only  the  magnitudes  of  both  fields  are  taken  into  account. 

In  computing  the  RI  noise  from  the  powerlines,  a  quasi-static  condition 
is  assumed.  This  is  valid  if  the  wavelength  of  the  frequency  under 
consideration  is  large  compared  to  powerline  parameters  (height  and 
phase  spacing)  and  to  the  distance  between  the  powerlines  and  the 
measuring  points.  Thus,  at  any  point  too  far  from  the  lines,  the 
validity  of  this  assumption  is  questionable.  However,  for  complete¬ 
ness,  low  ERP  and  laige  distances  are  being  considered  in  this  report. 
Fig.  2.13  to  2.24  have  indicated  that  for  a  range  of  low  ERP  and  the 
line  designs  considered,  if  the  location  of  the  powerline  is  close  to 
an  NDB,  the  critical  distance  is  small  at  any  receiver  altitudes.  Thus 
it  can  be  concluded  that  if  a  powerline  is  close  to  an  NDB  transmitter, 
the  RI  noise  radiated  by  the  powerline  will  practically  cause  no  appre¬ 
ciable  effect  to  an  aircraft  flying  past  it. 
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Fig.  2.13a  Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =0.1  watt,  line  voltage 
345  kV,  f  =  200  kHz,  under  heavy  rain  condition. 
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ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =0.1  watt,  line  voltage  ■ 
345  kV,  f  =  500  kHz,  under  heavy  rain  condition. 
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NDB  transmitter  to  powerline;  ERP  =0.1  watt,  line  voltage 
500  kV,  f  =  200  kHz,  under  heavy  rain  condition. 
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ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =  0,1  watt,  line  voltage 
500  kV,  f  =  500  kllz,  under  heavy  rain  condition. 


ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =0.1  watt,  line  voltage 
765  kV,  f  =  200  kHz,  under  heavy  rain  condition. 
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ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =  0.1  watt,  line  voltage  ^ 
765  kV,  f  =  500  kHz,  under  heavy  rain  condition. 


Fig.  2.16a  Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =  0.1  watt,  line  voltage 
1100  kV,  f  =  200  kHz,  under  heavy  rain  condition. 


ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
^^DB  transmitter  to  powerline;  ERP  =  0.1  watt,  line  voltage  ^ 
1100  kV,  f  =  500  kHz,  under  heavy  rain  condition. 


lAla  Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
N'DB  transmitter  to  powerline;  ERP  =  0.5  watt,  line  voltage 
345  kV,  f  =  200  kliz,  under  heavy  rain  condition. 
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.17a  Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
N'DB  transmitter  to  powerline;  ERP  =  0.5  watt,  line  voltage 
345  kV,  f  =  200  kHz,  under  heavy  rain  condition. 


Fig.  2.17b  Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  pouerline;  ERP  =0.5  watt,  line  voltage  ■ 
345  kV,  f  =  500  kHz,  under  heavy  rain  condition. 
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Fig.  2.18a  Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =  0.5  watt,  line  voltage 
500  kV,  f  =  200  kUz,  under  heavy  rain  condition. 


ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =  0.5  watt,  line  voltage  ^ 
500  kV,  f  =  500  kHz, under  heavy  rain  condition. 
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Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  pov.erline;  ERP  =  0.5  watt,  line  voltage  ■ 
765  kV,  f  =  200  kHz,  under  heavy  rain  condition. 


Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =0.5  watt,  line  voltage 
765  kV ,  f  =  500  k)Iz,  under  heavy  rain  condition. 


Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =  0.5  watt,  line  voltage  ■ 
1100  kV,  f  =  200  kHz,  under  heavy  rain  condition. 


ratio  as  a  function  of  aircraft  altitudes  and  distance  froi 
NDB  transmitter  to  powerline;  ERP  =  0.5  watt,  line  voltage 
1100  kV,  f  =  500  kHz,  under  heavy  rain  condition. 


NDB  transmitter  to  pouerline;  ERP  =  1.0  watt,  line  volta; 
345  kV,  f  =  200  kHz,  under  heavy  rain  condition. 


Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =1.0  watt,  line  voltage 
345  kV,  f  =  500  kllz,  under  heavy  rain  condition. 
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Fig.  2.23a  Critical  distance  f  i  oi.,  .;ircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  fiincti;^n  of  iircr.ift  altitudes  and  distance  from 
NDB  transmirti  r  to  po-.t-rl  i ne;  i  KP  =1.0  watt,  line  voltage 
7fi5  kV,  f  -  'ur)  Ml-,  o  !>  r  lu-avy  rain  condition. 


Chapter  III 

I’lUiDICTION  or  RADIO  INTliRFORCNCE  NOISE  PROM  DC  POWERLINES 


3.1  Introduction 

D i  rcct-current  powcrlines  can  be  clas.sified  into  three  types, 
Theyarc;  1)  the  monopolar  link,  2)  the  bipolar  link  and  3)  the 
homopolar  link. 

1)  The  monopolar  link  has  one  conductor,  usually  of  negative  polarity, 
and  ground  or  sea  return. 

2}  riie  liipolar  link  has  two  conductors,  one  positive  and  the  other 

nci;.'i  t  i  VC', 

3)  The  homopolar  link  has  two  or  more  conductors  all  having  the  same 
polarity,  usually  negative,  and  always  operates  with  a  ground 
return . 

RI  noise  from  DC  power  lines,  like  its  counterpart  in  AC  powcrlines,  is 
also  caused  by  corona  discharges  on  the  line  conductors.  Studies  and 
report  ^ have  consistently  shown  that  the  major  source  of  RI 
noise  is  from  the  positive  conductor  and  that  contribution  from  the 
negative  conductor  to  the  total  RI  noise  is  practically  negligible. 
Based  on  the  conclusion,  the  prediction  method  outlined  here  for  the 
RI  noise  from  the  DC  powerlines  will  only  consider  the  bipolar  link. 

As  opposed  to  RI  noise  from  AC  powerlines,  RI  noise  from  DC 
powcrlines  decreased  under  heavy  rain  conditions.  The  above  mentioned 
studies  and  reports  have  concluded  that  under  foul  and  wet  weather 
conditions,  the  RI  noise  level  decreased  considerably.  However, 
tlic  RI  noise  level  increases  under  fair  weather  conditions. 

Another  result  of  the  investigations  is  that  there  is  a  large 
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12 

effect  of  wind  velocity  on  the  level  of  RI  noise.  The  CPRl  test 
concludes  that  the  RI  noise  levels  arc  increased  by  wind,  with  the 
greatest  influence  being  when  the  direction  of  air  flow  is  from  nega¬ 
tive  to  positive  conductors.  A  wind  velocity  of  15  meter/second  (33.3 
mph)  will  cause  an  increament  of  about  5  dB/pV/m.  Another  report^"^ 
concludes  that  the  attenuation  of  the  RI  noise  level  laterally  from  the 
positive  conductor  is  not  appreciably  affected  by  conductor  configura¬ 
tion  or  pole  spacing. 

3.2  Computation  of  RI  Noise  from  DC  Powerlines 

To  date,  tlicrc  arc  OTily  two  empirical  formula  developed  for  the  calcu¬ 
lation  of  RI  noise  from  PC  powcrlines.  One  has  been  developed  by  EPRI 

14 

and  the  other  by  Reiner  and  Gehrig  .  Since  EPRI  is  the  authority  on 
this  subject  matter,  we  will  use  their  equation  for  the  prediction  of 
RI  noise  from  DC  powcrlines. 
riic  RI  noise  level  can  ho  calculated  by 

17  g  2 

I.  =  214  log -  278  (log  -  )  +40  log  a  (3-1) 

^o  ^o 

where:  i;  is  the  RT  noise  level  in  dB/  V/m  at  834  kllz  and  tlic  location 

of  the  observer  is  at  30.5  m  from  the  positive  conductor, 

g  is  the  maximum  conductor  surface  gradient  in  kV/cm, 

®max 

g^  is  the  critical  gradient  in  kV/cm, 
a  is  the  radius  of  subconductor  in  cm. 

This  is  essentially  the  maximum  value  of  the  RI  noise  level  as  the 
equation  is  formulated  for  fair  weather  condition.  For  completeness, 
correction  factors  for  frequency  and  radial  distance  arc  added  into  the 
above  cqua: ion.  Therefore,  more  generally 
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g  2 

n  =  21A  log  -  278 (log  +  40  log  a  + 

^o  ^o 

OT  ,  >IA  1  30.5  ,, 

27  log  —  +  40  log  (3-2) 

where  f  is  the  frequency  in  kilohertz, 

d  is  t!ic  radial  distance  from  the  positive  conductor  in  meters. 

liPRI  found  g^  to  he  approximately  14  kV/cm  for  the  conductors  tested. 

I'his  value  will  be  used  as  a  constant.  For  the  value  of  g  ,  we  will 

*’max 

use  the  equation  developed  by  Mangoldt^^.The  gradient  factor  for  the 
bipolar  DC  powerlines  is  given  by 


c  = 


1  ^■(  (N  -  a)  a/R  ) 


N  a  In 


211 


(3-3) 


R  ^  1 

eq  J  S 


where  (1  is  the  gradient  factor  in  kV/cm  per  kV  to  ground, 

K  is  the  radius  of  an  equivalent  bundle  conductor,  where 


eq 


l<  •  (Na  U"-')'''* 
eq  ^ 


(3-4) 

N  is  the  number  of  suheonductors  in  the  bundle, 
a  is  the  radius  of  subconductor  in  cm, 

R  is  the  radius  of  the  circle  on  which  the  centres  of  the  sub- 
conductors  lie  with  R  =  (B/2)/sin(ti/N]^where  B  is  the  distance 
between  tlic  adjacent  subconductors, 

II  is  the  average  height  of  conductors,  defined  as  the  minimum 
height  above  ground  plus  1/3  of  the  sag, 

.S  is  the  pole  to  pole  spacing. 

The  maximum  bundle  surface  gradient  is  thus  obtained  by  the  equation 

g  -  C  •  V  (3-.S) 

i'tax  ' 


IS 


! 


where :G  is  the  gradient  factor  in  kV/cm  per  kV  to  ground, 

V  is  the  line  voltage  to  ground  in  kV. 

EPRI  has  also  made  an  observation  that  the  RI  noise  levels  attenuate 
ns  the  square  of  the  rndinl  distance  until  it  roaches  an  inflection 
point  beyond  which  the  attenuation  is  proportional  to  the  radial 
distance.  The  inflection  point  can  be  approximated  by  the  equation 


H  ^ 
^  -  2^ 


(3-6) 


where  d  is  the  radial  distance  from  the  positive  conductor,  and  is 
the  wavelength  of  the  noise  frequency.  Beyond  this  point,  equation 
(.1.2)  can  bo  rewritten  as 


g  g  2 

V.  =  214  log  -  278  (log  )  +  40  log  a  ♦ 

®o  ®o 

,  834  ,  30.5 

27  log  j—  *  20  log  -g— 


(3-7) 
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For  illustration,  some  actual  line  designs  will  be  considered.  Table 
3.1  details  the  line  voltages  and  parameters.  Fig.  3.1  to  3.4  show 
the  characteristics  of  the  calculated  RI  noise  levels  for  each  of 
the  powerline  designs  listed  in  Table  3.1,  at  the  frequencies  of  200 
kHz  and  500  kHz  with  the  location  of  the  observer  at  various  altitudes. 


Line 

Config. 

Line  Voltage 
(±  kV) 

H(ft) 

s(ft) 

a  (cm) 

N 

n 

400 

63.00 

34.50 

3.048 

1 

#2 

450 

40.00 

44.00 

2.032 

2 

113 

600 

49.87 

36.75 

1.525 

4 

U 

750 

54.40 

45.00 

2.032 

4 

Table  3.1  DC  line  voltages  and  parameters  considered 


where: H  is  the  average  height  of  the  conductors, 

S  is  the  spacing  between  the  conductors, 
a  is  the  radius  of  the  subconductors, 

N  is  the  number  of  subconductors  in  a  bundle. 

Appendix  E  details  program  DCRI  listing  for  the  computation  of 
the  RI  noise  level  radiating  from  the  DC  powerlines. 

3.S  Prediction  of  Critical  Distance  Where  the  Ratio  of  Desired 
Signal/Undesired  Noise  is  15  dB 

Similar  to  Section  2.4,  the  desired  signal  is  the  signal  from  the 
NDB  transmitter  (covered  in  Section  2.3),  and  the  undesired  signal 
is  the  RI  noise  from  the  DC  powerlines  dealt  with  in  Section  3.2. 
Consider  the  same  scenario  shown  in  Fig.  2.12  in  Section  2.4. 

Fig.  3.1  to  .‘?.4  have  shown  that  the  levels  of  the  RI  noise  radiated 
from  the  DC  powerlines  are  lower  than  that  radiated  from  the  AC  power- 
lines.  Therefore,  only  line  ±400  kV  which  radiates  highest  RI  noise 
level  at  higher  altitudes  is  considered  here.  Fig.  3.5  and  3.6 
indicate  the  critical  distance  for  this  line  design  at  various  receiver 
altitudes  with  the  F.RP  of  the  NDB  being  0.05  watt  and  0.1  watt  respec¬ 
tively. 

3.4  Conclusion 

As  in  the  case  of  the  AC  powerlines,  a  quasi-static  condition  is  assumed 
and  the  computation  of  the  ratio  of  desired  signal/undesired  noise 
involves  the  magnitudes  of  the  fields  only.  For  illustration,  Fig.  3.5 
and  3.6  show  the  critical  distances  for  an  NDB  transmitter  with  an  ERP 
O.Ofi  watt  and  0.1  watt  respectively.  These  values  of  ERP  are  too  low 
for  any  practical  purposes.  It  can  be  predicted  that  at  a  higher  range 
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of  ERP  values,  the  critical  distance  will  reduce  almost  to  zero.  Thus 
if  a  bipolar  DC  powcrline  is  located  close  to  an  NDB  transmitter,  the 
RI  noise  radiated  by  the  poworlino  should  cause  no  serious  effect  to 
an  aircraft  flying  past  it. 
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at  '30  kHz  under  fair  weather  condition. 


Calculated  RI  noise  profile  for  a  bipolar  ;400  kV  DC  powerline 
at  500  kHz  under  fair  weather  condition. 


Fig.  3.2a  Calculated  RI  noise  profile  for  a  bipolar  ±450  kv  DC  powerline 
at  200  kHz  under  fair  weather  condition. 
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Fig.  3.3b  Calculated  R1  noise  profile  for  a  bipolar  ±600  kV  DC  powerline 
at  500  kHz  under  fair  weather  condition. 


at  200  kHz  under  fair  weather  condition. 


Calculated  RI  noise  profile  for  a  bipolar  t750  kV  DC  powerline 
at  500  kHz  under  fair  weather  condition. 


ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =  0.05  watt,  line  voltage 
±400  kV,  f  =  200  kHz,  under  fair  weather  condition. 


Fig.  3.5b  Critical  distance  from  aircraft  to  powerline  for  15  dB  s/n 
ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =  0.05  watt,  line  voltage 
±400  kV,  f  =  500  kHz,  under  fair  weather  condition. 
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ratio  as  a  function  of  aircraft  altitudes  and  distance  from 
NDB  transmitter  to  powerline;  ERP  =0.1  watt,  line  voltage  ■ 
±400  kV,  f  =  200  kHz,  under  fair  weather  condition. 
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Chapter  IV 

RliRADIATION  PREDICTION  FROM  POWER-LINE  STRUCTURES 

4.1  Introduction 

In  order  to  provide  predictions  of  the  possible  effects  on  NDB 
reception  of  signals  reradiated  by  large  conducting  objects,  such  as 
power  transmission  lines,  accurately  and  economically  a  large  scale 
series  of  measurements  is  not  feasible.  In  order  to  provide  a  large 
enough  data  base  for  meaningful  conclusions  a  great  many  airborne 
measurements  of  field  strength  would  need  to  be  made.  Furthermore, 
actual  sites  which  would  provide  suitable  situations  regarding  the 
relative  locations  of  an  NDB  and  large  power  line  structure  would  have 
to  be  located. 

Fortunately,  there  has  been  considerable  progress  made  recently 
in  tbe  application  of  computers  to  electromagnetic  problems.  More 
specifically,  techniques  have  been  developed  which  enable  the  calcula¬ 
tion  of  induced  currents  and  the  resulting  reradiated  fields  on  wire 
structures  caused  by  incident  electromagnetic  energy  provided  that  the 
wire  structures  are  small  in  comparison  to  the  wavelength  of  the 
incident  field.  Since  NDB  frequencies  are  so  low,  even  very  large 
mechanical  structures,  such  as  sections  of  power  transmission  lines, 
are  small  in  comparison  to  the  NDB  wavelengths. 

Since  these  computer  techniques  have  not  been  previously  applied 


to  the  situations  we  are  here  concerned  with,  a  validation  by  comparison 
with  measurement  seemed  appropriate.  Measuring  the  reradiated  fields 
from  an  actual  power  line,  however,  would  not  be  possible  due  to  the 
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very  complex  geometry  of  the  structure.  Instead,  a  radio  tower  has  been 
used  for  the  validation.  Since  the  source  of  reradiation  from  a 
straight  vertical  tower  is  localized,  the  reradiated  fields  can  be 
measured  by  using  a  loop  receiving  antenna  and  orienting  the  loop  so 
as  to  null  the  signal  coming  directly  from  the  NDB.  If  the  geometry 
is  chosen  suitably,  the  loop  may  simultaneously  be  oriented  so  as  to, 
at  least  approximately,  respond  maximally  to  the  reradiated  fields  com¬ 
ing  from  the  tower.  A  suitable  tower  is  one  owned  by  radio  station 
WRrn  in  Columbus,  Ohio.  Tt  is  a  straight  tower  558  feet  high.  The 
NDB  transmitters  in  the  vicinity  are  CMH  and  DKG.  The  former  is  oper¬ 
ating  at  391  Kliz  with  an  effective  radiated  power  of  9.683  watts  and 
the  latter  at  .348  KHz  with  an  HRP  of  0.6  watt.  Fig.  4.1  is  a  map 
showing  the  locations  of  the  tower  and  the  NDB  transmitters.  The  CMII 
antenna  is  situated  at  about  11.72  NM  S-23®-E  of  the  tower  and  that  of 
DKG  is  at  about  5.86  NM  S-27®-W.  A  sketch  map  is  illustrated  in  Fig. 
4.2. 

4.2  Measured  Results 

Measurements  were  taken  along  the  access  road  leading  to  the  tower. 
Equipment  used  included  Interference  Analyzer  Model  EMC- 25,  Remote 
Vertical  Antenna  Model  RVR-25  and  Remote  Loop  Antenna  Model  ALR-25. 
Measurements  were  made  at  five  different  places  along  the  access  road 
as  indicated  by  Fig.  4.3.  In  each  case  the  vertical  or  loop  antenna 
was  mounted  on  a  five- foot  high  tripod. 

Results  obtained  were  the  total  electric  field  strength  as  mea¬ 
sured  liy  the  vertical  antenna  and  the  magnetic  field  strength  as  mea¬ 
sured  liy  tlie  loop  antenna.  When  using  the  loop  antenna  two  measurements 


Radio  tower  (558') 


1  inch  =  100' 


Fig.  4.3  Sketch  of  the  locations  of  five  measuring  points 
with  respect  to  the  radio  tower  and  the  NDB 
transmitters. 


were  made.  The  loop  was  rotated  in  a  vertical  plane,  and  the  maximum 
and  minimum  values  of  the  magnetic  field  were  measured.  The  maximum 
value  was  considered  to  be  primarily  due  to  direct  radiation  from  the 
NDB,  and  is  denoted  the  "direct"  magnetic  field.  The  total  electric 
and  direct  magnetic  measured  fields  are  listed  in  Tables  4.3  and  4.4. 

4.3  Calculated  Results 

The  calculated  results  are  obtained  by  using  a  computer  program 
based  on  the  work  of  Richmond^^.  The  program  makes  use  of  a  mathemati 
cal  model  of  the  NOR  transmitting  antenna  and  the  radio  tower  to 
compute  the  electric-field  strength  for  any  location.  The  ground  is 
assumed  flat  and  perfectly  conducting,  as  the  major  interest  is  the 
rcrad i a t ed  si gna 1 . 

r.ssentially,  the  purpose  of  this  preliminary  work  is  to  chock 
the  accuracy  and  reliability  of  using  this  particular  computer  program 
to  predict  signal  reradiation.  The  test  involves  comparison  of  the 
measured  results  and  the  calculated  results  that  are  obtained  by  using 
the  mathematical  models  that  simulate  the  same  physical  situation. 

The  coordinates  of  the  five  measuring  locations  with  respect  to 
each  transmitting  antenna  are  shown  in  Tables  4.1  and  4.2.  In  each 
case  the  x-axis  is  the  imaginary  line  that  passes  through  the  NDB 
antenna  and  the  reflecting  tower  with  the  former  being  the  origin. 

The  electric  field  strength  of  the  direct  signal  from  the  NDB  trans¬ 
mitters,  the  reradiated  signal  from  the  tower  and  the  total  electric- 
field  strength  are  calculated  along  the  access  road  with  respect  to 
CMII  .and  DKCi.  These  are  illustrated  in  Fig.  4.4  and  4.3.  In  addition 
Fig.  4.(>  shows  the  calculated  behavior  of  the  total  electric-field 


Location 

X- coordinate 
(feet) 

y-coordinate 

(feet) 

z-coordinate 

(feet) 

II 

1 

71.065.0 

-71.7 

5.00 

# 

2 

71,020.0 

-42.3 

5.00 

« 

3 

70,980.0 

-16.1 

5.00 

# 

4 

70,895.0 

39.6 

5.00 

5 

70,815.0 

92.0 

5.00 

Table 

4.1.  Coordinates  with  respect  to  CMH  transmitter 
The  x-axis  originates  at  the  NDB  and  passes 
through  the  reflecting  radio  tower. 

Location 

x-coordinate 

(feet) 

y-coordinate 

(feet) 

z-coordinate 

(feet) 

II  1 

35,557.0 

-168.0 

5.00 

II  2 

35,509.5 

-185.6 

5.00 

»  3 

35,462.0 

-203.2 

5.00 

#  4 

35,367.0 

-238.3 

5.00 

II  5 

35,269.5 

-274.4 

5.00 

Tabic  4.2.  Coordinates  with  respect  to  DKG  transmitter. 

The  x-axis  originates  at  the  NDB  and  passes  through 
the  reflecting  radio  tower. 
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along  tlic  x-axis  for  +300  feet  from  the  tower  for  the  CMH  transmitter, 
fig.  4.7  illustrates  the  similar  characteristics  for  the  DKG  trans¬ 
mitting  antenna.  In  both  cases  the  curves  are  smooth  as  expected  and 
the  total  electric-field  strength  levels  increase  at  points  closer  to 
the  tower  due  to  the  fact  that  the  reradiated  signal  strength  is  higher 
in  this  region.  The  peak  value  of  the  tower  electric-field  strength 
that  coincides  with  the  location  of  the  tower  should  be  ignored  as  this 
value  i.s  not  valid.  The  computer  model  cannot  handle  points  inside 
a  conductor. 

4.4  Comparison  of  the  Measured  and  Calculated  Results  for  Direct  and 

Total  Fields 

Moa.sured  results  obtained  by  using  the  vertical  antenna  arc  com¬ 
pared  with  the  calculated  total  electric-field  strength.  Fig.  4.4a 
and  4.5a  show  the  levels  of  the  total  electric-field  strongtli  for  the 
CMII  and  DKII  transmitting  antenna  respectively,  and  comparison  between 
the  measured  and  calculated  results  can  be  made  directly.  For  conven¬ 
ience,  Tables  4.3  and  4.4  provide  the  comparison  for  the  five  measuring 
locations. 

In  the  calculated  results,  the  electric-field  strengtli  of  the 
direct  signal  from  the  transmitting  antenna  could  be  converted  to  a 
magnetic-field  strength  by  using  the  equation  p  =  E/H,  when  n  is  the 
intrinsic  impedance  of  the  medium  and  could  be  assumed  to  lie  approxi¬ 
mately  120  .  The  measured  levels  of  the  magnetic-field  strength  (with 
loop  oriented  for  maximum  signal)  for  the  CMH  and  DKG  transmitters 
arc  also  shown  in  Fig.4.4h  and  4 . 5b  rc.spcctively .  Comparisons  arc 
given  IS  well  in  Tables  4.5  and  4.6. 


Calculated  and  measured  values  of  electric  field  strength 
along  the  access  road  with  respect  to  CMH  transmitter. 


Fig.  4.4b  Calculated  values  of  radiated  magnetic  field  strength  and 
measured  maximum  magnetic  field  strength  along  the  access 
road  with  respect  with  CMH  transmitter. 


Fig.  4.5b  Calculated  values  of  radiated  magnetic  field  strength  and 
measured  maximum  magnetic  field  strength  along  the  access 
road  with  respect  with  DKG  transmitter. 


Calculated  values  of  the  total  electric  field  strength  with 
respect  to  DKG  transmitter  along  the  x-axis,  originating  at 
the  NDB  antenna  and  passing  through  the  tower. 
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Local  ion 

Measured  Results 
(dB/yV/m) 

Calculated  Results 
(dB/yV/m) 

Difference 

(dB/yV/m) 

#  1 

69.6 

63.3 

6.3 

#  2 

68.6 

63.0 

5.6 

#  3 

68.6 

62.8 

5.8 

«  4 

68.6 

62.6 

6.0 

#  5 

68.6 

62.5 

6.1 

Table 

4.3.  Total  electric-field  strength  comparison  for 

CMIl  transmitter.  Measurements  made  using 
vertical  monopolc  antenna 

Location 

Measured  Results 
(dB/yV/m) 

Calculated  Results 
(dB/yV/m) 

Difference 

(dB/yV/m) 

It  1 

59.8 

57.30 

2.. 50 

#  2 

59.8 

57.05 

2.75 

#  3 

59,8 

56.90 

2.90 

»  4 

59.3 

56.70 

2.60 

It  5 

58.8 

56.65 

2.15 

Tabic  4.4,  Total  electric-field  strength  comparison 
for  UKG  transmitter.  Measurements  made 
using  vertical  monopole  antenna. 


Location 

Measured  Results 
(dB/yA/m) 

Calculated  Results 
(dB/yA/m) 

Difference 

(dB/yA/m) 

#  1 

11.9 

11.165 

0.735 

«  2 

11.4 

11.170 

0.230 

»  3 

11,9 

11.175 

0.725 

#  4 

11.9 

11.186 

0.714 

II  5 

11.9 

11.196 

0.704 

Tabic  4.5,  Direct  magnetic-field  strength  comparison 
for  CMH  transmitter.  Measured  results 
are  for  loop  oriented  for  maximum  response. 


Location 

Measured  Results 
(dB/yA/m) 

Calculated  Results 
(dB/yA/m) 

Difference 

(dB/yA/m) 

H  1 

-2.4 

5.20 

-7.60 

It  2 

-1.9 

5.20 

-7.10 

H  3 

-1.4 

5.22 

-6.62 

II  A 

-1.4 

5.24 

-6.64 

II  5 

-2.4 

5.25 

-7.65 

Table  4.6.  Direct  magnetic-field  strength  comparison  for 

DKG  transmitter.  Measured  results  are  for  loop 
oriented  for  maximum  response. 


radiation)  along  the  access  road  with  respect  to  DKG 
transmitter. 
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4.5  Near- Zone  Magnetic-Pield  Strength  Comparison  of  Calculated  and 

Measured  Rcradiated  Fields 

Wien  the  loop  antenna  used  to  measure  magnetic  field  strength 
at  the  access  road  measurement  points  was  rotated  to  the  minimum  signal 
position  for  the  DKG  signal,  the  direct  signal  from  the  NDB  was  approxi¬ 
mately  in  the  loop  null.  At  all  the  five  measuring  locations,  the 
alignment  of  the  loop  was  simultaneously  such  that  strong  coupling 
existed  with  the  currents  induced  on  the  tower. 

The  values  of  the  magnetic  field  measured  then  were  due  primarily 
to  the  rcradiated  signal  from  the  tower,  other  interferring  signals, 
and  noise.  Unfortunately,  this  situation  did  not  happen  for  the  CMII 
transmitting  antenna.  Thus  only  DKG  signals  are  compared  in  this  sec¬ 
tion. 

For  near-zone  rcradiated  magnetic-field  strength  calculation, 
a  further  separate  computation,  also  based  on  the  work  of  Richmond^^^^ 
was  necessary.  The  calculated  results  arc  shown  In  Fig.  4.8  and 
comparison  with  the  measured  values  is  listed  in  Table  4.7. 


Location 

Measured  Results 
(dn/uA/m) 

Calculated  Re.sults 
(dB/uA/m) 

Difference 

(dB/uA/m) 

II  1 

-9.4 

-19.4 

10.0 

II  2 

-11.4 

-20.. 

R.9 

II  ^ 

-1.'^.4 

-21.2 

7.R 

II  4 

-16.4 

-2.%.1 

6.7 

II  5 

-20.4 

00 

1 

4.8 

Table  4.7  Near-zone  rcradiated  magnetic-field  comparison 
for  the  DKG  transmitter.  Measured  results  arc 
for  the  loop  oriented  for  minimum  response. 
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The  base  of  the  tower  is  insulated  from  ground,  and  therefore  the  calcu¬ 
lated  values  in  Table  4.7  were  made  for  this  condition.  The  measured 
results  are  consistently  higher  than  these  calculated  results  indicating 
that  more  current  was  induced  in  the  actual  tower  than  predicted  by  the 
computer  model.  The  actual  tower  is  fed  at  the  base  through  a  matching 
device  which  provides  a  current  path  from  the  tower  base  to  ground.  How¬ 
ever,  this  impedance  could  not  be  conveniently  measured  since  the  radio 
station  broadcast  24  hours  per  day,  nor  did  the  station  engineer  have  this 
information.  However,  the  calculations  were  repeated  with  the  tower  base 
shorted  to  ground,  and  both  open  and  short  circuit  results  arc  shown  in 
Fig. 4. 8.  The  measured  results  for  the  actual  tower  lie  nearly  midway 
between  the  two  calculated  curves.  Since  the  problem  is  linear,  one  would 
expect  that  if  the  correct  tower- to-ground  impedance  were  input  to  the 
model,  the  discrepancy  would  necessarily  be  less  than  indicated  in  Fig. 

4.8,  which  is  on  the  order  of  10  dB.  Also,  since  transmission  line  towers 
are  well-grounded,  this  ambiguity  will  not  exist  in  the  intended  applica¬ 
tion,  and  the  conclusion  is  that  the  computer  modeling  program  is  valid 
for  predicting  reradiated  signal  levels  at  NDB  frequencies. 

4.6  Prediction  of  the  rcradiatcd/dircct  signal  ratio  for  various 
powcrline  tower  models. 

At  this  point,  we  arc  able  to  investigate  the  levels  of  reradiated 
signals  from  various  powcrline  tower  models  by  using  the  computer 
modeling  program.  For  this  purpose,  three  different  tower  configura¬ 
tions  are  being  considered.  Their  detailed  dimensions  are  shown  in 
Fig.  4.9.  Each  tower,  in  turn,  is  placed  at  1  NM  mile  away  from  the 
NHR  transmitter  along  the  imaginary  x-axis,  which  is  passing  through 
the  transmitter  and  the  powerline  tower  with  the  former  being  the 


Fig.  4.10  A  general  layout  of  a  powerline  tower  model  with  respect 
to  an  NDB  transmitter. 
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origin,  as  illustrated  in  Fig.  4.10. 

The  reradiated/direct  signal  ratio  is  computed  along  the  x-axis 
within  ±200  feet  from  the  tower  and  the  receiver  altitudes  are  at 
500,  1000  and  1500  feet  above  ground  respectively.  The  value  of  BRP 
of  the  NDB  transmitter  is  immaterial  for  this  purpose  as  only  the 
ratio  of  the  reradiated/direct  signal  is  being  investigated.  Any  value 
of  l=RP  of  the  NDB  transmitter  will  be  cancelled  out  in  the  computation 
for  the  reradiated/direct  signal  ratio.  Similarly,  the  distance 
separating  the  NDB  transmitter  and  the  powerline  tower  does  not  affect 
the  ratio  greatly  either.  In  fact,  as  the  separating  distance  in¬ 
creases,  the  reradiated/direct  signal  ratio  decreases. 

The  levels  of  the  reradiated/direct  signal  ratio  for  the  three 
different  tower  models  at  the  frequencies  of  200  and  500  KHz  arc  shown 
in  I'ig.  4.11,  The  plots  consistently  show  that  the  rerad iated/dircct 
signal  ratio  is  much  lower  than  the  critical  value  of  -15  dB.  It 
appears  that  the  different  powerline  tower  models  do  not  affect  the 
rcradiated/direct  signal  ratio  a  great  deal,  as  indicated  by  Fig.  4.11. 
Based  on  this  observation,  a  simplified  powerline  tower  model,  as  shown 
in  Fig.  4.12,  can  be  conveniently  used  for  further  investigation.  Fig. 
4.13  confirms  that  the  reradiatcd/dircct  signal  ratio  of  this  simpli¬ 
fied  powerl Lne  tower  model  is  comparable  to  those  of  the  three  different 
tower  models. 

Using  this  simple  tower  model,  a  general  layout  as  shown  in  Fig. 
4.14,  can  be  considered.  Three  tower  models  are  placed  1000  feet  apart 
and  their  alignment  is  perpendicular  to  the  x-axis  at  1  NM  away  from 
the  NDB  transmitter.  The  resulting  rerad iated/direct  signal  ratio  is 
shown  in  Fig.  4.15.  The  plots  indicate  that  the  reradiated/direct 


Fig.  4.11a  Reradiated/direct  signal  ratio  for  double  wire  steel  tower 
model  at  200  KHz  vs  distance  from  the  NDB.  The  tower  is 
at  6076  feet  (i  NJ4)  and  the  flight  path  is  directly  over 
the  tower. 


model  at  500  KHz  vs  distance  from  the  NDB.  The  toi 
at  6076  feet  (1  Wl)  and  the  flight  path  is  djtectly 
the  tower. 


Fig.  4.11c  Reradiated/direct  signal  ratio  for  single  steel  tower  model 
with  cross  arms  at  200  KHz  "s  distance  from  the  NDB.  The 
tower  is  at  6076  feet  (1  MM) and  the  flight  path  is  directly 
over  the  lower. 
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Fig.  4. lie  Reradiated/direct  signal  ratio  profile  for  double  steel  tower 
model  with  cross  arms  at  200  KHz  vs  distance  from  the  NDB. 

The  tower  is  at  6076  feet  (1  N>;)  and  the  flight  path  is 
directly  over  the  tower. 


The  tower  is  at  6076  feet  (1  W4)  and  the  flight  path  is 
directly  over  the  tower. 
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.13a  Reradiated/direct  signal  ratio  profile  for  the  simple  tower 
at  200  KHz  vs  distance  from  the  NDB.  The  tower  is  at  6076 
feet  (I  NH)and  the  flight  path  is  directly  over  the  tower. 
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Fig,  4.15a  Rerad iated/direct  signal  ratio  profile  for  the  three  tower 

model  layout  at  200  KHz  vs  distance  from  the  NOB,  The  tower 
is  6076  feet  (1  t€1)and  the  flight  path  is  directly  over 
closest  tower. 
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signal  ratios  are  well  below  the  critical  value  of  -IS  dB. 

Further,  possible  effects  of  the  horizontal  conductors  are  also 
investigated.  The  conductors  are  modeled  as  a  single  eight- inch 
diameter  smooth  conductor,  spanning  2000  feet  across  the  three  simple 
towers,  as  shown  in  Fig.  4.16.  The  resulting  reradiated/direct  signal 
ratio  is  plotted  in  Fig.  4.17.  As  anticipated,  the  reradiated/direct 
signal  ratios  of  the  horizontal  conductor  model  are  much  lower  compared 
to  those  computed  for  the  vertical  tower  models.  The  reason  is  that 
the  vertically-polarized  radiated  signal  from  the  NDB  transmitter 
reacts  minimally  with  the  horizontal  conductors. 

4.7  Conclusion 

The  results  of  Section  4,6  have  consistently  shown  that  the  reradiated/ 
direct  signal  ratio  is  well  below  the  critical  value  of  -15  dB  even  for 
flight  paths  only  500  feet  above  ground.  Thus,  it  can  be  concluded 
that  the  reradiated  signals  from  powerline  structures,  line  conductors 
included,  should  not  interfere  with  normal  ADF  operation. 


Fig.  4.16  Model  of  horizontal  line  conductor  spanning  across  the 
three  towers  with  respect  to  an  NDB  transmitter. 


The  tower  is  at  6076  feet  (1  Kl) and  the  flight  path  i 
directly  over  the  closest  tower. 
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Fig.  4.17b  Reradiated/direct  signal  ratio  profile  for  the  horizontal 
line  conductor  model  at  500  KHz  vs  distance  from  the  NDB. 

The  tower  is  6076  feet  (1  Nil) and  the  flight  path  is  directly 
over  the  closest  tower. 
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CONCLUSIONS 


The  objective  of  this  report  is  to  assess  the  possibilities  of 
harmful  interference  with  ADF  operation  being  caused  by  electric  power 
transmission  lines.  The  two  modes  of  interference  considered  are  corona 
radio  noise  generated  actively  by  the  transmission  lines,  and  passive 
reradiation  of  the  radio  signal  from  the  NDB  by  the  transmission  line 
structures. 

In  order  to  achieve  this  objective,  computer  programs  were  devel¬ 
oped  which  can  predict  the  critical  distance  between  the  aircraft  con¬ 
taining  the  ADF  receiver  and  the  power  transmission  line  where  the 
ratio  of  desired  signal  to  either  transmission  line  noise  or  reradiated 
signal  is  l.S  dB.  These  predictions  can  be  made  as  a  function  of  the 
distance  separating  the  NDB  transmitter  and  the  powcrline,  the  elevation 
of  the  flying  aircraft  (and  hence  the  ADF  receiver),  the  ERP  of  the  NDB 
transmitter,  the  geometry  and  voltage  of  the  powerline,  and  the  frequen¬ 
cy.  For  the  case  of  generated  radio  noise,  the  permittivity  and  conduc¬ 
tivity  of  the  earth  can  also  be  varied;  however,  for  the  passive 
rcradiation  case  the  earth  is  assumed  to  be  perfectly  conducting. 

Regarding  active  generation  of  radio  noise,  AC  lines  generate 
considerably  more  radio  noise  than  do  DC  lines.  However,  even  for 
the  worst  case  considered  in  this  study,  a  765  kV  AC  powerline  located 
20  nautical  miles  from  a  200  kHz  NDB,  the  results  show  that  for  an 
aircraft  flying  at  only  a  500  feet  altitude,  the  aircraft  if  on  a  path 
through  NDB  location  and  perpendicular  to  the  powerline,  would  have  to 
approach  to  loss  than  4000  feet  from  the  powerline  before  the  ADF 
signal  to  noise  ratio  dropped  below  15  dB.  The  re.sult  is  for  a  weather 
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condition  of  heavy  rain,  which  causes  an  AC  line  to  generate  the  highest 
level  of  radio  noise.  If  the  N'DB  is  located  more  closely  to  the  power¬ 
line,  the  radio  noise  has  relatively  less  effect  on  ADF  performance 
since  the  desired  NDB  signal  level  is  increased  while  the  radio  noise 
level  remains  unchanged. 

Moving  the  NDB  closer  to  the  power  transmission  structures  does, 
however,  increase  the  likelihood  of  ADF  interference  due  to  passive 
reradiation  of  the  NDB  signal  by  the  powerline  structures.  The  possi¬ 
bility  of  such  interference  was  estimated  using  a  computer  simulation 
model  based  on  the  method  of  moments.  However,  even  for  the  worst  case 
considered  which  had  a  500  kHz  NDB  located  one  nautical  mile  from  a  single 
steel  tower  with  the  aircraft  flying  over  the  tower  at  an  altitude  of  500 
feet,  the  reradiated  signal  was  calculated  to  be  approximately  38  dB  Lzlow 
the  level  of  the  directly  radiated  (desired)  signal.  The  fundamental 
reason  for  this  is  that  the  wavelength  of  the  NDB  signals,  even  at  500 
kHz,  is  large  compared  to  the  tower  dimensions. 

Thus  the  conclusion  is  that  locating  an  NDB  near  a  high  voltage 
transmission  line  (up  to  within  1  nautical  mile)  would  probably  not  impair 
the  function  of  the  NDB  due  to  either  corona  noise  generation  or  passive 
reradiation  from  the  line.  However,  as  shown  by  the  worst-case  results 
previously  mentioned,  there  is  a  possibility  of  interference  for  a  situ¬ 
ation  which  involves  a  low  power  beacon  located  at  a  relatively  large  dis¬ 
tance  from  a  potentially  interfering  power  line.  And,  it  should  be 
pointed  out  that  other  possible  interference  mechanisms,  such  as  power¬ 
line  carrier  radiation,  have  not  been  considered  in  this  study. 
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APPENDIX  A 


Derivation  for  equation  of  field  factor 
for  point  P  anywhere  in  space 
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=  60  •  1  • - —  (A-9) 

x2  +  h2 

The  direction  of  the  electric  field  is  vertical,  perpend iculc'  to  the 
magnetic  field  as  shown  in  Fig.  A-1. 

Normally,  when  dealing  with  magnetic  field,  it  is  adequate  to  consider 
the  line  conductor  in  free  space  without  any  image.  However,  if  the 
earth  is  assumed  perfectly  conducting,  the  image  will  act  to  double  the 
field  at  ground  level  due  to  symmetry.  In  all  the  sections  dealing 


with  powcrline  conductors  in  this  report,  it  is  assumed  that  the  earth 
is  perfectly  conducting  and  thus  fields  due  to  image  conductors  are 
also  being  considered. 

For  any  point  in  space,  as  shown  in  Fig.  A-2,  the  magnetic  field  due  to 
corona  generated  currents  in  the  powerline  conductor  can  be  determined. 
The  horizontal  component  of  the  magnetic  field  is  given  by 


H  =  H.cosB  -  Hpcos  a 
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Similarly,  the  vertical  component  of  the  magnetic  field  is  given  by 

(A-14) 


H  =  H„sin  a  -  H. sin  3 
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itiuh,  Liio  t'csi.ltani.  nuRnclic  field  is  given  by 
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whose  direction  is  as  shown  in  Fig.  A- 2. 
The  corresponding  electric  field  is  then 
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Simplifying  equation  fA-21)  yields 
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The  direction  of  tne  electric  field  is  perpendicular  to  the  magnetic 
field  as  shown  in  .'ig.  A- 2.  liquation  (A-22)  can  be  written  as 

P.  =  1  •  F  (A-23) 


whore  F  is  termed  the  field  factor. 

in  applying  modal  analysis  for  the  computation  of  the  RI  noise  radiated 
from  the  3-phasc  AC  powerlines,  the  resultant  RI  noise  for  each  phase  is 
related  to  tlio  modes  by  the  equation 


where  m  denotes  the  mode  number  and  n  =  1,2  or  3.  This  operation  is 
detailed  in  section  2.2. 
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1  line  conductor 


Vie.  \-l  Illustration  of  the  magnetic  and  electric  fields 
®  at  ground  level  due  to  current  I  in  a  line  conductor 

;ibovc  a  perfectly  conducting  ground. 


r. 

1 


rig.  A-2.  Illustration  of  the  magnetic  and  electric  fields  at  any 
point  in  space  due  to  current  I  in  a  line  conductor 
above  a  perfectly  conducting  ground. 
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APPliNUIX  B 


Program  CHARGE  listing  for  computing  the  values 
of  the  total  charge  on  each  conductor  bundle 
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')ATA  XL/*D*  ,*T*  .•  A*  .•N»,»C***E»#»I*  t*F»,*T»,*  )•/ 

^ATA  Pl.EPSICn/l* 14 I59.a, R54F-12/ 
iM  A  FPSI  ,S!r,MA,E»P.C/!0.aO,O.Ol00.0.50,1.0E8/ 

HE  16,54) 

hnPM4r(«i«,TlO,«N*,Tl6,*KV* , t 26, •hT • * T16, *0* , T46, •« AO* , TS5, • R« AO* 
R  t  A-)  15,1)  N,  V,H,r),  R  AO,  A 
FrmMAT(tio,5Fio.i » 

<PfrF  16,1)  N,V,H,0,AAD,A 
REAO  15,2)  I  on  ) ,  !«  I,  1) 

<:nPM4V<?^ \S,1) 

WHME  16,52) 

i-nfi'*Arf*;)S5X,*VALUFS  OF  CHARGE  ON  EACH  CONDUCTOR*) 

WR(1<^  (6,2)  (0(ll,(«t,T> 

1.1  ntrumNE  THF  MAXIMUM  BUNDLE  GRADIENT  FOR  EACH  PHASE 
writs:  ((,,S0I 

FnhMAT|*0»,5X, ‘VALUES  OF  MAXIMUM  BUNOtE  GRADIENT*) 

•;N  -  N  -  1 
on  10  la  1 , 1 

PFO  «  (ELOAriN)*BAO*IAPANN)l*A( l/N) 

^  AVI  -  CArtSIOI! 5 )/(  ^.OOPPf*FPSILOAftOAT|N»ARAOAl.OOE1) 

F’JAVI  n  =  EAVf  *<  I  ,no*l  FLOAT!  NNIPR  AO/A)  ) 

HPI  ^F  <f.,1)  I  .ewaxi  r  J 

EnF«‘MTUOX,*EMAX  *•  ,F10.  3, 2X  ,  •  KV/CM  •  ) 

LfiNT  f  NUE 

lij  OETFR«|NE  THF  EXCITATION  FUNCTION  ITAUI  OF  EACH  PHASE 
(iNHER  HFAVV  RAIN  CONDITION 
WRfTF  (6,511 

fnPMAT(*0*,5X, •VALUES  OF  EXCITATION  FUNCTION*) 

TFIN. EO. I »  KN  -  T 
ir(N.F0.2)  KN  »  2 
IFIS.GE.1)  KN  «  0 

U(A  =  2.00BRA0 

on  20  I  *  1,1 

TAUOR  a  76,00-(550.00/EMAX(lllB3B«00«AtOG10IDIA/3.90l«^FLOAT1KNI 


acroooio 

ACR00020 
ACROOOn 
ACR00040 
AC« 00050 
ACR00060 
ACROOOTO 
ACRO^lflO 
ACB00090 
ACROO) 00 
ACROOl 10 
AC»0f)l20 
ACROOl in 
ACROOl 40 
ACROO) SO 
ACRO)  160 
Af ROOl TO 
ACROOl RO 
Afo  00  190 
Afi»00?00 
ACR0O2I0 
ArRn02?0 
ACRno210 
)ACR00240 
ACR00250 
ACR00760 
ACR00270 
ACR002RO 
ACR00290 
ACROOIOO 
ACROOIID 
ACR00320 
ACROOIIO 
ACR00140 
ACR00150 
ACR00160 
ACR00170 
aCROOIRO 
ACROOIBO 
ACR00400 
ACR00410 
ACR00420 
ACR00410 
ACR00440 
AfR 00450 
ACR00460 
ACR00670 
ACROQARO 
ACR 00490 
ACR00500 
ACROOSIO 
ACROr»520 
ACR00530 
ACR 00540 
ACR00550 


138 


FORTRAN  A  OHIO  UNIVERSITY  Fl6C  fROHAGNf  T  IC  C.OH*»A  T  AR  |  L  I  F  IF$  LAB 


rt  .<  :  TAUn‘^/?0.00 

T  4IJ1  !  )  -  I  10,00**RA  T» 

^^rrf  ifTAuin 

4  Fn‘MATnOX,"TaU  «•  ,F  10.^,2X.  •HfCROAMP/SORF  N*| 

?0  C’^^riNUF 

c 

C  TO  OETERNIN6  THE  'GFQHETRIC  MATRIX*  OF  THE  PHASES 

^0 ! TP  I4f S1I 

Si  Pni<HAT|*0*.SX, •VALUES  OF  GEOMETRIC  MATRIX  IC>*> 

'TO  ^0  I  -  U3 
o  ri  n  *  H 
I  (  n  =  -H 
no  TO  J  *  1,3 
I  r-  I  I  .eo,  J  I  GO  TO  31 
nw( ! , j»  =  0 

I  E  M  ,E0.  I  .  ANO.  J*EQ.  ^  I  OwM.Jl  *  r)Y2*00 
•..F').3.4Nn,J.fO.  11  OWil.JI  *  f'*?.00 
f)  1  i,J)  a  SOOTKMTI!)  HKin**?  ♦  0WII*J»*R2> 

G!(,J)  a  ALQOIOHll, JI/OMf ftJll 
GO  TO  30 

51  G(I,J1  *  ALOGI?.OORHTniAlOO*00/A| 

30  continue 

v4«irF  lA.SI  ((Gn,J)  ,J«l,3lf  I«1,3) 

*5  1  nPMATn<SX,FlO,3ll 

c 

c  if)  iiETERMINE  THE  •MODAL  TRANSFMATfON  MATRIX*  OF  THE  PHASES 

C  MtEO  TO  COMPUTE  TmE  VALUES  OF  LAMflOA  WHERE  LAMROA  ARF  THE  CURIC 

C  ROOVS  FOR  THE  DETERMINANT  OF  •GEOMETRIC  MATRIX*  EQUAL  ZERO 

C  ♦  Gll-LAM  G12  G13  • 

C  *  GM  G22-LAM  G23  ♦  *  0*0  lAI 

C  •  GU  G32  G33-LaM  R 

C 

C  OBTAIN  the  CUflIC  EQUATION  ANO  SOLVE  FOR  LAMBDA  ILAM) 

C  the  cubic.  EQUATION  IS  OF  THE  FORM 

"  I  AM{cUBEI^Al*LaMIS0>*A2»LAM*A3  *  0*0 

C  TO  SOLVE  EOR  LAMBDA 

WRITE  IN, SSI 

SS  Et!R«AT  (  •  0»  ,sx,  *  VAIU^^  >  OF  LAMBDA*! 

A!  -  (GliWl  »  r,l2,cM  •  013,31) 

A?  ^  (  G(  I ,  I  I  *0  I  ?,  21  ♦  r,(  I ,  u  *01  T,  31  »  GI?*2IPGf3f  31)  - 
/  (  II  ♦  G(  3,?)Ar,I?,3»  ♦  Gn.3)RGI3f  I)) 

A3  ^  (Gf  I,  n*G(  3,21 ‘GI  2,  31  ♦  G I  U  ?)  ♦G<  ?*l  IPGI  3,3 )  ♦ 

?  Gri,3)««G(3,l)*G(2,2l»  -  I G  f  1 ,  U  *01 PGI 3,  3)  ♦ 

3  Cl  1,?)RGI3,IIMGI2,3)  ♦  G I  1 , 3)  RGf  ?« 1 1*01  3, 2M 

00  ^  <3,OOMA2  -•  AloAll/9,0a 

R  =  I9.00PAI*A2  -  2T.OO*A3  -  2.00*1 AIR*3» 1/54.00 
TMEVA  a  ARCOSIR/SORTI-Oon 
DO  40  I  *  1,3 

IftMi  I  )  a  2.00*SQRT(-00)*CnSMTHETA/3.00»*’IFlOATf  II-1.00)*I20.00* 
2  PI/180.00I-AI/3.00 

HOITE  iNf?)  IfLAMU  ) 

>  FORMAH  lOX,  *L  AM  (•,11,*)  *«,E10.3I 

•♦0  continue 


ACROOS40 
ACROOS70 
ACROOSAO 
ACRO0S9O 
ALR00600 
ACROOMO 
ArR00620 
ACR00430 
ArRD0640 
AC *006 SO 
ACR00440 
ACR00470 
ACROOftRO 
AC.  R  00490 
ACROO7O0 
ACROOTtO 
Af.R007?0 
ACR00730 
ACR00740 
ACROOTSO 
AC» 00740 
ACP0O770 
APR007AO 
ACR00790 
ACROOROO 
ACR008I0 
ACR00320 
ACR00R30 
ACR00B40 
ACROORSO 
ACROOB60 
ACR00R70 
ACR00B90 
ACROOB90 
ACR00900 
ACP009I0 
ACR00920 
ArR00930 
ArR00940 
ACP009S0 
ACR00940 
AC.ROQ970 
ACR009BO 
ACR00990 
ACR01009 
ACROtOlO 
ACR01020 
ACR0I030 
ACROIOAD 
ACR01050 
ACROIOAO 
ACROI070 
ACR01080 
ACR0I090 
ACROltOO 


Fite:  ACRI 


FORTRAN  A 


OHfO  UNiVFRStTY  ELFCTROHAGNETIC  CO'IP A T AA II !  T I E S  LAB( 


c 

AC901 110 

r 

TO  »R))*NGE  THE  VALUES  OF  LAMnOAS  IN  ASCFNDI>4G  OROER 

AC801 120 

WRITE  (IS,S6I 

AC801 130 

56 

FORWATI'G'.SK, 'VALUES  OF  LAHROAS  IN  ASCENDING  DROER'I 

AC80I 160 

1  =  J 

AC80tl50 

NN  -  I-l 

AC801160 

DO  AO  K  >  1,NN 

ACROIITO 

JJ  *  I-A 

AC8  0I  180 

on  TO  L  -  l,JJ 

Ar.«01  190 

IFILAHILI.LE.LAMIL'LII  go  to  to 

AC801200 

TFHR  -  LANILI 

ftC80l?10 

lAMILI  «  LANIL'll 

AC«017?0 

LAWILFll  ■  TEMP 

AC801230 

ro 

CDNIINUF 

Arft01260 

60 

CONTINUE 

AC801250 

WRITE  <A.«I  II.LAMI 1 l,|3l,3| 

AC801760 

8 

FORNATI lOX.'LAN  ('.ll.'l  »',F|0-AI 

AC901270 

C 

AC80I 280 

C 

SURSH  LUTING  EACH  OF  THESE  VALUES  OF  LAMBDAS  INTO  EQUATION  lAI 

AC801290 

c 

WILL  GIVE  US  THE  COLUHN  OF  THE  'MODAL  TBANSFORMATION  MATRIX'  INI 

AC80I 300 

WRITE  IS.STI 

Ar.Bomo 

37 

FORMATI'O'.SX, 'VALUES  OF  MODAL  TRANSFORMATION  MATRIX  |M|'| 

Af:8ai32o 

C 

FOR  FIRST  COLUMN 

AC891330 

DO  AO  1  >  1,3 

ACR0I34O 

80 

Gll.ll  «  Gll.ll  -  LAMfii 

AC801330 

Mil. II  <  1.00 

AC80t360 

Ml x,l 1  •  Nit, 11 

AC801370 

MI?, II  •  -IGI 1,1I*MI1,I|  ♦  GII,3I*H|3,III/CI 1,21 

ACR0I380 

OIVl  •  S0RTIm|I,||«*2  ♦  Nl?,l>**2  ♦  Mt3,'l*R2l 

AC801390 

DO  »l  I  •  1,3 

AC801400 

81 

M|t,l  1  >  M(l,ll/OIVl 

AC931410 

C 

ACR01420 

c 

FOR  SECDNO  COLUMN 

ACR01639 

on  998  1  •  1,3 

AC801640 

998 

Gll.ll  «  Gll.ll  *  LAMIll  -  LAM(2I 

AtROlASO 

Ml  1,71  «  -1.00 

Ar«0!660 

MU, 21  .  -Mll,?l 

ACROIATO 

M|?,?l  .  -IGI  t.l  l*MI  1,21  *  GUt3l*MI3,2||/GU,2l 

AC»01680 

0IV7  •  SQRTfl  l,2l**2  ♦  MI2,2I»»7  ♦  MI3,2I**2I 

AC801690 

on  99t  I  •  1,3 

ACR01800 

997 

Ml  I,T|  «  MII,?I/0IV2 

Ar«0l5l0 

C 

ACR01520 

c 

FOR  THIRD  COLUMN 

AC801830 

00  8?  1  •  1.3 

ACROIBAO 

82 

Gll.ll  «  Gll.ll  F  LAMI2I  -  LAMI3I 

AC801880 

Mil, 31  <  1.00 

AC8  01860 

"(3,31  •  "11.31 

AC801870 

"12,31  -  -(G(l,ll*M(l,3l  ♦  G(t,3|RM(3,3MTGIl,2) 

AC* 01 980 

niV3  •  SQRTIMI 1,3I**?  ♦  M(2,3|R»7  ♦  MI3,3I*«2I 

AC801990 

no  83  I  -  1,3 

AC801600 

83 

"11,31  •  Mn.3l/DIV3 

ACR01610 

C 

AC801620 

WRITE  IA,9I  ll"ll,JI,J«l,3l,(«l,3) 

AC801630 

9 

FnRMATI3ISX,Fin.3ll 

ACtt0l640 

C 

ACB0I690 

F  1 1  i  :  4  r  P  I 


OHIO  UNIV6PSITV  fLFCTRnHAGNFTIC  COMPAT AR II  I T I  6 S  lAB 


c 

S 

Af  POl  <S60 

c 

Tn  nPTERHiNf  IHF  Rnfiii  cnMpnufNts  <if  cnRn><«  ru««FN>  injfctio'is 

ACPni6T0 

c. 

NEFO  TO  FIND  THF  INVEASf  HATAICfS  Of  IGI  ANO  Ml 

ACR  OI6AO 

c 

ACR0I69D 

c 

TO  FINO  INVF9SF  MATRIX  Of  IGI  ICII 

ACRniTOO 

WH|Tr  (6,Se) 

AFROlTin 

SB 

FnRM'\TI«0**SX,»V4LOFS  OF  |NV€RSe  •MATRIX  Of  IC»*I 

ACROl  770 

no  1  «  l#3 

AC»01 730 

R4 

GM  .  n  s  G(  I.  1  1  *  LAHI  31 

ACROl 740 

X  =*  1 

AC^O! 760 

'lALL  iNVfRS  lK,G,M,;i| 

ACROl  760 

00  90  J  ■  1,3 

ACR0177f) 

on  00  1  *  1,3 

AC0017R0 

'iO 

r.  It  1 ,  J )  *  ;  111 «  J I 

ATroI 790 

werTF  (6,t?i  nGn[,j),j-u3i,t  =  i,3i 

ACROl ROO 

r 

ACROl  RIO 

C 

TO  FINO  INVMSe  MAIRIX  OF  |M|  |.F*  MU 

ACROl B70 

-‘M  r  t  f  6  ,  S9I 

ACROl  ft30 

>  i 

n  »  M  .  f  •  C  * ,  s  . ,  •  /44  iF  S  OF  1  NVfRSF  MATRIX  OF  M)  •  | 

ACR0I840 

K  »  ? 

ACROl PSO 

K.til  INV6RS  (K,r.,M«/t| 

ACROl P60 

nn  91  j  ^  l,  1 

ACROIR70 

on  91  1  •  1,3 

ACR  OIRRO 

M  (  (  ,  J  )  «  Z 1  (  1  ,  J  1 

ACROl  R90 

WPITF  <6,121  1(9111, 3), |*t,1> 

ACR01900 

\2 

FM«MAT<3(SX,FlO,ni 

AfROlRIO 

C 

ArR0l970 

C 

rn  noTAiN  thf  ppoouct  of  <gu*mti 

ACROI930 

UR] TF  (6,61) 

ACROIRAO 

FnRMAr(»0»,SX,*VA(.a«:S  OF  PROOUCT  OF  (MIIRIGIIM 

ACR0I9S0 

nn  9?  r  1,3 

ACR01960 

-■)0  9  •»  J  ^  1,3 

ACR0I970 

GM< I , JI  a  0,0 

ACR019R0 

;jn  92  <NOgX  »  l»  ] 

ACR01990 

-  GH(I,JI  4  Mf(l,|NOFX|PGMINOEX,J) 

ACR02000 

ON  t  (  NUH 

ACR07010 

f6fl2»  ((CMU  ,j|,j-»l,3),l  =  l,3l 

ACR02020 

c 

ACR02n30 

c 

ACR02040 

c. 

.  «SS«»S3S«S« 

ACR02050 

c. 

(  r.M'MJT  ION  FOP  PHiSC  ONE  ONLY 

ACR02060 

c 

.  i.  :  ^s.-:jw;sss=W'-«>-.s33.sxaS9«3as«3 

ACR02070 

c 

ACR070R0 

c 

TO  ORVAIN  THE  VALUES  OF  THF  CORONA  CURRENT  INJECTION  JC 

OH 

ACR02090 

r. 

PHaSE  ONE 

ACR02I00 

WRITE  (6,141 

ArR02ll0 

14 

FnRMAT(«0*,SX, 'VALUES  OF  CORONA  CURRENT  IN  EACH  MODE’ 1 

ACR021?0 

00  93  1  *  1,3 

ArR02l30 

93 

JCI  I,  l»  *  CN(f  ,1  IPTAUd) 

ACR02140 

WRITF  (6,991  (jri f • 

ACR02190 

99 

FORMAT! 15X,F10.3I 

ACR02160 

C 

ACR02I70 

r 

VO  OFTIftHINE  THE  MQOAL  CURRENTS  TN  THE  THREE  CONDUCTORS 

FOR 

ACR02190 

t 

PHftSF  ONE 

ACR02I90 

WR  n  6  (6, ISI 

ACR02200 

OHIO  (l*4fVfRS|TY  f  Lf  C  f  ROM  4GNF  f  I C  COMRA  TAB  IL  I  T  T  f  S  I  ABf 


FORMATI'O*  ,*5*#  • 'VALUES  OF  HOO^L  CURRENTS  IN  EACH  CONOUCTOPM 

f  •  ?00,00 

CAtL  CONST  <V,F,At<»HA| 
on  R4  I  •  l,B 

IC2lfl«l)  *  JCU  ,  I) /l?.00*S0«TT  ALRHAf  in  ♦ 

WRITE  I6,16)  (l.ic?in«n»f  =  l»3) 

FORM  ATI  IS  X,  •  lf.?l  I  *  t  f  U  =•  .F  10.  3  ,  ?X,  •  Ml  CRPAMP*  ) 

F  «  MOO.OO 

CALL  CONST  (V.F,AtPHAI 
no  9S  I  •  1,3 

icsiiiti)  «  jci  f«i)/(2.oo«softrf  ALPHAtnn 

WftlTF  16,131  I  I  ,  irSH  (•  tl  ,  l  =  t 

format  1 1 sx, • icB I  ( • , 1 1, • • n  =• .FI 0.3,2 X, • -ir voamp*  > 

TO  OFTERMINF  THF  PHASE  CURRENTS  |P|I,J) 

FOR  FREOUENCV  AT  200  KH? 

WRITE  f6,Tll 

FORMAT|*0*,SK, ‘VALUES  OF  MODAL  CURRENTS  IHICROA)  AT  203  KHM) 
K  •  I 

CALL  IPMASE  CM ,K , KK ,KKX, IC21 . If SI , IC22, IC52, I C2 3, I C S3 , I p I 
00  96  I  «  1,3 
00  96  J  »  1,3 
|P2tl  l,JI  -  IPCf ,J) 

WRITE  16,121  inP2lf!,JI,J=l, 31, 1^1,31 

FOR  FREOUENCV  AT  SOO  KHZ 
WRITE  16,721 

FORMATI*OSSX, ‘VALUES  OF  MODAL  CURRENTS  IMICROA)  AT  SOO  KHZM 
K  «  2 

CALL  IRHASE  C  M  ,K  ,  KK  ,  KKK,  I C  2 1  ♦  ICS| ,  I C22,  tCS2,  IC  23,  1C  S3 , 1 1»  I 
00  97  I  *  1,3 
on  97  J  *  I  ,3 
fRSlI  UJ)  «  IPC  If  Jl 

WRITE (6, 12)  Cl IRStl I, J), 4*1,3), 1*1,3) 


COMOtJTATION  EOR  PHASE  TWO  ONLY 


WRITF  C6,16» 

00  72  f  *  1,3 

Jf C I , 1 )  *  GMC I ,2) •TAUC7I 

WRITE  16,99)  f Jf f I, I )•{*! ,3) 

TO  DETERMINE  MODAL  CURRENTS  IN  THE  THRFF  CONOUTTORS 
WRITE  16, IS) 

F  «  200.00 

CALL  CONST  <V,F, ALPHA) 
on  23  T  ■  1,3 

IC2?t1,l)  ■  jcn,l)/f  7.00*SQRTCAlPMAfin) 

WRITE  16,26)  ( I, IC22C I, n ,|«),3) 

FORMAT! ISX, • If 22  «• ,F 10. 3, 7X, • MI CROAMP • I 


ACPrj27lO 
ArP0?22O 
ArPO2?30 
Ar.p  02260 
ArfiO?2SO 
ACR02260 
A''»O?2?0 
ACR022A0 
Afoa2290 
Ar»02300 
Ar.R02310 
Afo 02320 
ACR02340 
ACRO73A0 
ACR023SU 
ArR02360 
ACR02370 
ACR J21A0 
tr9n?^Q0 
Art-  02600 
ACRn26  10 
afp  02620 
ACR02630 
ATP  0  2660 
ArRO26S0 
ACR02660 
ACR02670 
ACR026R0 
AfR 02690 

Acon?soo 
ArRO?AlO 
ACR02620 
ATROZSTO 
iCRors60 
ACRO2SS0 
4FP07S60 
ACR02ST0 
ACR02SRO 
ArPD?SO0 
ACR  02600 
AfRO?6lO 
ArRa?b?o 
Ar»026^0 

4r.a02660 

ACRO26S0 

ACR02660 

ACR02670 

ACR02630 

ACR02690 

ACR07700 

ACR02T10 

ACR02720 

ACRD2730 

ArR02760 

ACR07TS0 


u  u 


■ 


14 


file:  *CRI 


FORTRAN  A 


OHIO  UNIVfl»*;ltV  ELFCTRONAGNEIIC  COH*A  I  AR  1 L  I  T  If  S  lAR' 


F  »  500.00 

r.Alf  CONSTtV.F.AlPHAl 

on  ?5  I  »  1,5 

?S  Tr5?(I,ll  -  JC( I, ll/l 2.00«S0«rC»tPHAIII »l 
WHTF  I  I.  IC52<  I.  M 

?6  FORMATT 15*, • IC5?  « • ,FI 0. 5. 2X, •MICROAMP' ) 

TO  OFTFRNINE  THE  PHASE  CURRENTS 
FOR  FRCOUENCY  at  200  KH2 
WRITE  (6,711 
RK  *  1 

CALL  (PHASE  (H,K,KK,KKK,  IC7l,(C5H,IC22,  IC52,  IC23,  IC53,  l»» 
on  32  (  -  1,3 

on  32  J  *  1,3 

32  IP2?(I.J1  •  IP(I,JI 

WRITF  (6,121  ( ( !P22( I , Jt , J*t,3t ,(=l.3» 

C 

c  FOP  FREOUENCY  at  500  KHZ 

HOITE  (6,721 
KK  =  2 

call  (PHASE  (H,K,KK,KKK,(C2l,ir5(,IC22.IC52.(C23,IC53,|P| 
no  36  I  >  1,3 

on  16  j  »  1,3 

36  |P5?(I,JI  ■=  (P(I,JI 

WRITE  (6,121  ( ( (P52( I,JI ,J=1.3l.(»l.3l 

C 

c 

^  3  as  ««  a  sasassa  aaaa  aaaaaaaaaaaaaaa 

c  tO‘^»UTfcT10N  THRtf  ONtY 

C 

WR(TP  f6«  141 
on  !0l  I  ■  1*3 

!f)l  JCdfll  •  GM(f,3l#TAU(3l 

WRire  fft*99J  (JC(  I* 

c 

c  TO  DfTPRMfNF  THF  M004L  CURRENTS  IN  THF  THRfF  CONDUCTORS 

WRITE  16, IS) 

F  s  ?00.00 

rail  CONST  <V,F,ALOHftl 
00  10?  I  a  1,3 

10?  ic?3fi,i)  -  jci r, n/( ?,oo*soRrf AtRHAf im 
WRITE  16,101)  ri,IC?3(I,U,(«l,3t 

103  FORNATI  ISX,*  IC?3  «  •  ,F  I  0.  3, 2X,  •  MKROA^P  •  ) 

c 

F  a  SOO.OO 

CACt  CONST  IV,F,ALPHA> 

00  104  la  1,3 

104  1C53n,ll  -  JCU,!I/I2,00*SORTCALPMA|IIII 
WRITE  (6,105)  ( I «  IC53I r*l )*I«I,3I 

105  F0RNATI15X,MC53  »•  ,F|  0*  3,  ?X, ’Nf  CROANP  •  ) 

C 

C  TO  OETFRHINF  THF  THF  PHASE  CURRENT 

C  FOR  FRFOUENr.Y  AT  ?00  KH? 

WRITE  (6,711 


ACRO?760 
A'-»o?7ro 
ACRO?  750 
Ar»0?79(J 
ACRO?  ROO 
ACQ0?R10 

ACO0?53O 
«^CR  0?540 
ArP0?850 
AC»0?860 
Ar«0>870 
Ar 

ACRO?%90 
ACR  0?900 
ACRO’OlO 
A^R0?9'>0 
ACi>0?<J30 
Ar!5r.>Q40 

Ar;/0>960 
ACR0?970 
ACR  07930 
ACR0?990 
ACROIOOO 
Af'oaioio 
ACR030?0 
ACR03030 
ACR'^1040 
ACR01050 
ACRa3060 
ArR03070 
Ar.R030R0 
ACR03090 
ACR03100 
ATROII  10 
arRon?o 
ATROll 30 
arR03!40 
ACR03150 
ACR03160 
ACR031T0 
ACR03180 
Ar.R03I90 
BCR03200 
ACR03?10 
ACR03??0 
ACR03?30 
ACR03?40 
ACR03?50 
ACR03?60 
ACR01?70 
ACR03?R0 
AC»01?90 
ACR03300 


uuu 
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FILE:  »CRI  fO»T«*N  L 


nmn  univeositv  elec  runmcNEitc  roMPitso  ii  i  r  |t  ^  L»»r 


106 

c 

c 


LOT 

c 

c 

c 

i; 

C 

C 

C 

c 

c 

c 

c 

c 

TOO 


62 

C 


C 


c 


i 


KKK  >  I 

CALL  IPHASt  IM,K,KK  ,KKK,  !C?l.  ICSl.IC??,  ICS2,  ICJ»  ICS^i  |P» 

DO  106  I  «  l«1 

DO  106  J  «  1«3 

|P21(  t»J)  «  t(*(  I  f  Jl 

wpfrr  (6,121  (UP2i( 

FOR  FRfOUFNCV  AT  500  KHl 

URITF  (6,72) 

KKK  *  2 

CALL  IPHASF  (M.K.KK.KKK, IC7l,fC5I .IC226 rC52, fC25, /C53,(PI 
on  107  I  •  1*3 

on  (07  J  *  1,3 

|P53( I ,JI  «  IPI I • Jl 

WRITE  16,12)  l(IP53(l«J),J»l,3),t^l,3) 

GO  TO  700 


PART  (A)  MILL  CO*<PUTE  THE  Rl  NOISE  LEVEL  FROM  THE  AC  TRANSmISION 
LINES  ALONG  THE  LATERAL  OISTANf.E  WITH  THE  RECEIVER  ALTITUDE  AT 
CROUNO  level. 

PART  IB)  WILL  COMPUTE  THE  CRITICAL  OfSTANCES  WHERE  THE  kAT(0  fiE 
OFSIRFO  SlftNAL  (FROM  NOR ) /UNOFS IREO  NOISE  IRI)  IS  15  DR  WITH  ThE 
RECEIVER  ALTITUDES  AT  500*,|000*  AnO  1500*. 

s  ssaaatsasssssa*  ssass33sasssasss&ssss«s9ss  ses  s  xxss  sssassssssss 


PART  (A) 

CONTINUE 

WITH  KNOWN  NO04C  CURRENTS  CALCUtATEO  ABOVE  THE  RI  NOISE  AT 
GROUND  LEVEL  COULD  RE  OETERNIKEO  NOW 
FOR  FREOUENCV  AT  200  KHZ  AND  500  KHZ 
WRITE  16*62) 

FORNATI»OSTlT,  'YLOC*  ,T30f  •XLnc*tT44,*EN?00»  ,  TSRt'ENSOO*  I 
lALTI  «  0 

VLOC  *  FLOAT!  IALT  MPO.IOAR 
XF  *  500.00 
XI  *  0.0 
NP  »  100 

XO  «  (XF-Xn/FLOAT(NP) 

XP  •  XI 

00  77  I  «  l,NP 
XP  s  XP  »  XO 
J  »  I^NP 

XLOC  *  XPP0.30N9 
K*l 

CALL  NOISE  (K,XtOC,YtOC,M,O.IP2l*IP5l*IP22,IP52,!P?3,IP53,EE) 

FA2  *  FE 

K»? 

CALL  NOISE  CK,XLOC,YLOr,H,D,IP21*lP5l,|P22*lP52,IP23, IP53.EE) 
EA5  *  FE 

VIII  -  FA2 


A  r; » ri  i  ^ )  0 
Ar»,033?0 
Ara033)0 
aCR  /33*0 
ACROTTSO 
ArR03360 
ACR03  370 
ACR  033«0 
4r.Rn33QD 
ACB  03400 
aCRf>34)0 
ArP03420 
aCP' 0  3430 
arp  )3‘*40 
AC  k O34S0 
AC»0346C 
ACR034‘?0 
ACR034B0 
ACR03490 
irft03500 
A  r  ft  0  3  I  0 
ArPD3S.>0 
ACR0353' 

Af.R03540 

ACR035S0 
ir.R03560 
ACP035Tf) 
Aro03580 
4rp^1^90 
ArR03600 
ArR036lO 
ACR  03620 
47003530 
arRO‘»,6<.o 
AC" 03650 
arR03660 
a'‘P  03670 
ACR036S0 
4CR036P0 
AC«01  700 
Ar.R037l0 
aCR03720 
ACR03730 
Arpn3740 
ACR03750 
ACR03760 
aCR03770 
AC»037RO 
ACRO3790 
ACR03800 
ACRO38I0 
ACR03870 
ACRn3R30 
ACR03R40 
ACRO305O 
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r^iiF:  Arm  fortran  a 


OHIO  UNIVERSITY  EtFC  IROMACNE  T  I C  COHP  A  T  AB  I  L  I  T  I  (  »■  I  API 


XCI 1  »  XP 

ATP  0^400 

vij)  •«  Fas 

acrotato 

XCJI  *  XP 

ACROXP^O 

c 

ACPOTSOO 

witrrf  (6,76)  Miri,  xri),y(i>«ir(ji 

ACROT900 

FORMATClOXt l t8,lFlS.3l 

ACP0T910 

77 

C0NT1NUF 

AC«0X920 

CALL  MPLOT(X,V,NP,?,XLtVL) 

ar.P0>9S0 

GO  TO  900 

ACP039<.0 

C 

AfPOTOSO 

C 

ACP0T960 

c 

PART  (8* 

ACP0T970 

800 

CONTI NUF 

ACP  OT9flO 

c 

LFT  MFIGHT  of  ANTFNNA  8F  40*.  X*l  IF  FREOUFNCY  IS  ?00  XM7  AND 

ACROT9-;fO 

c 

Ka?  IF  frequency  is  SOO  KM/ 

Arp04ooo 

WRITF  (6,SS0» 

ACB0401 0 

sso 

FO«MAT(*O*,Tll,«00(NM|»,T24,»0l0)ST4O,«O(S00)*,TS6,»'J(ir»00I', 

Ar.Pi>4020 

2  T72,«O(tS00IM 

arc 

c 

ACPOSO^O 

K  »  2 

ACP040S0 

HI  »  40,00t0.1049 

ACRO4060 

XF  >»  20,00 

Ar«04070 

XI  «  0,00 

ArP04040 

NP  a  20 

ACR04O90 

OX  ■*  (XF-Xf  l/FLOATINP) 

ACP  04  100 

XP  a  X( 

Ar#04l 10 

on  400  I  *i,NP 

ACP04120 

XP  a  XP40X 

ifPOA  1  TO 

J  *  1  ^NP 

&rp04 140 

1  *  J4NP 

4C.P041S0 

N  *  t ♦NP 

ACPOA 1^0 

f 

Ar«04l 70 

GO  TO  SOS 

ACP041  a:» 

c 

FOR  RECEIVER  ALTITgoE  AT  GROUND  iFVFt  (H/aO.OM 

aroo4i90 

on  soo  INC  *  1,10000 

ACPO4200 

CALL  Sir.NAL  (  1 NC  •  K,  XP  ,Hl  ,  H2,  P|  •EPSI  ,  SI  GNA  ,  FRP,C  tf  F  S  1 

4CC:042  10 

SIGI  «  FFS 

4rff042.’O 

Yi-^r  »  H2 

Af.P042T0 

xinr.  *  float(inciao,io49 

arBn4240 

CALL  NHfSF  ( K , XLOC , YLOC • H, 0, (P21 1 1 PS  1 , rP22t 1 P52. 1 P23f f PS  1 , F E 1 

ArP04250 

RINI  *  EE 

ACP  04760 

OFITAI  a  StGWRINl 

Ar»0427C 

(FfoFlTAUGE^lS.OO)  CO  TO  600 

ar.F042P0 

soo 

Cr^NTINUE 

ACP  04?90 

600 

OISTI  •  FLOAT! INCI 

ACR04T00 

c 

ACP  04’.  10 

80S 

CONTINUE 

ACR04320 

c 

FOR  PFCEIVER  ALTITUDE  AT  500*  !M2«500-0*I 

Ar.R04TTO 

M2  a  S00.00P0.S049 

ACR04T40 

DO  SOI  INC  »  1, 10000*10 

ACR043Sn 

CAll  SIGNAL  1 INCfK,XP,Hl,H?,P|,EPS|,S(GNA,ERP,C,EFS) 

ACR04S6O 

SfG2  a  FFS 

ACP04T70 

nnr,  »  h2 

ArPOSTAO 

XtnC  *  FLnATIINC»P0,3049 

ArP04T90 

CALI  NOISF  (K,XLnr,YLaCtH*D*IP21*|P5t,IP?2,IPS2*IP2StlPST,EEI 

ArP044oo 
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FILE: 

*r»I  FORTRAN  A  OHIO  UNIVFRSItT  FtEC TROHAGNFt 1C  COHPA TAR  11 1 1 1 FS  Lt 

«IN2  ■  EE 

ACR04410 

DELTA?  -  SIG?<-R|N? 

' ACX04420 

IFIOFl TA?,GE. 1^.001  cn  T060t 

ACP04430 

SOI 

COST  1 HUE 

ACR04440 

601 

01  ST?  «  FLOATt IHC ) 

ACR044S0 

C 

ArR04460 

c 

FOR  RECEIVER  ALTITUOF  AT  1000*  |H?*1000*I 

Ar,R044TO 

H?  «  I000.00R0.1049 

ArR044flO 

DO  SO?  INC  «  !• 10000*10 

ACR 04490 

CALL  SIGNAL  ( 1 NC * K, XP* HI  * H?* PI tEFSl * SIGMA.f RP*C .FFS > 

ArRf>4S00 

SIGS  »  EFS 

ACR04S10 

VLOC  »  H2 

ACR04S20 

RLOC  3  FLOAT! INC) *0. )049 

aCR04SB0 

CALL  NOISE  |K.XLOC*YLOC*M*0*IP?l*IFSI,fP22*fPS2*!P2B*|PS1.EE) 

ACR04S40 

RINI  »  EE 

ACR04SS0 

0ELTA1  *  SICT>RtN3 

ACR04560 

IFCOELrAS*GE.lS.OO)  GO  TO  60? 

ACR04570 

so? 

CONTINUE 

ACR04S^0 

602 

OISTS  »  FLOATIINCI 

ACR04S90 

C 

&r.R04  600 

c 

FOR  RECEIVER  ALTITUOF  AT  1500«  fH2*lS00*l 

ACR04610 

H?  s  lS00*00R0«S04O 

AC«046?0 

00  SOB  INC  «  1*10000*10 

ACR04610 

CALL  SIGNAL  IINC* K, XP*Hl .H?* »I tEPSt • SIGHA* ERP*C • EFS 1 

ACR04640 

SIG4  *  FFS 

ArR046S0 

YLOC  «  M? 

ACR04b60 

XLOC  ■  FLOAT! INC) *0. 3049 

ACROAATO 

CALL  NOISE  <K*XLQC*YLQC*H*0*IP?UtPSl»IP229  lPS?*|P?3*fPS3*EEI 

ACR04640 

RtN4  «  EE 

ACR04690 

DELTAS  «  S1C4-RIH4 

ACR04TOO 

IF(0€LTA4.G6*15*00)  GO  TO  60S 

ACR047I0 

SOS 

CONTINUE 

ACR04r20 

60S 

0IST4  •  FLOATIINC) 

ACR04T30 

C 

ACR0474n 

YIII  «  XP 

ArR047X0 

X(|)  *  OISTl 

ACR04760 

rij)  *  XP 

Ar.!^  04770 

XIJI  »  OIST? 

ACRO47R0 

Y(L)  *  XP 

ACP04790 

XIL)  -  0IST3 

AfR  04900 

YIN)  *  XP 

AfPOASlO 

X(N)  s  0IST4 

ACR04RZ0 

c 

arR049xo 

c 

&CP04R40 

WRITE  (6,4011  Y(l 1,XIJ1*X(LI,X(NF*SIG?»R1N?,SIG3*RINS 

ACROASSO 

401 

FORNATISFS.l I 

ACR04R60 

400 

continue 

ACR048TO 

c 

CALL  NPLOT  I  X* Y,NP *  4  *  XL  * YL 1 

ACR04R80 

WRITE  16*1991  K,ERP 

ACP04  890 

194 

FnRNAr(SX*TS»FI0«1l 

ACR04900 

f 

ACR049I0 

GO  TO  2000 

ACR049?o 

C 

ACR04930 

t 

ACR04940 

900 

CONTINUE 

ACP049SO 

n  -o 
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ftlf-l  iCRl  fnUTHAN  A  OHtn  UNIVewSfTy  f-LfCtR0H4GNETIC  r  OMP  ATAft  !  L  !  T  16  S  LARf 


C 

C  THIS  PART  Will  fOMPlJTE  fHF  PI  NOISE  AT  ANV  POINT  OF 

C  OftSpPVATION, 

K  *  ? 

KK  «  2 

on  1000  I  ■  2,400,KK 
IFI I.GE.lOni  KK  >  2S 
Xinc  ■  PLOATU 

c 

H2  s  00(1.00P0.1049 
VLOC  •  H2 

CALI  NOISE  (K,XLOCfYLOC*H,0*|P2l*fP5l.|P22,fP52,IP23*IPSS.EP» 
EO  -  EE 
C 

H2  *  S(}0.00P0.3049 
Vine  »  H2 

CALL  NOISE  (K,XLOC#VLOCtH,0,IP2l.|P5l.IP22.IP5?,IP?3.IP^3.EEI 
ESOO  ■  EE 
C 

H2  «  t000.00A0«304<» 

Vine  *  H2 

CAIL  NOISF  IK,XL0C,rL0C,H,0tIP2l,IP5l#IP22«!PS?*|P23.|PS3*F6l 
6  inoo  «  EE 
C 

YLOC  *  1500. OOAO, 3049 

CALL  NOISE  IK,XLnCfVLOC,H,0,IP2l,IP5l»IP22f|P52*IP23,tP53,EE) 
EISOO  •  EE 

r 

WRITE  (4,1200)  I .EO.ESOOtElOOOtEtSOO 
1200  FORMAT  I5X,|5,4(F10.3)) 

1000  CONTINUE 

2000  continue 

STOP 

ENO 

C 

SURROUTINF  fNVFRS  IK, G, 9,21) 

01  mens  ION  G(30t30),2(30,30)f ZI(30«30),ZCf30f 30),ZA(  30,30) 

real  moo, 30) 

no  9  J  »  1,3 
00  9  I  ■  1,3 

IF(k.EQ.I)  2(1, J)  >  GdfJ) 

1FIK.EQ.2)  2( I,  J)  «  HI  I, J) 

CONTINUE 


C 


C 


TO  OBTAIN  COFACTOR  MATRIX  (2Cf 
Zrn,l)  ■  ♦(2(2,2)42f3,3)-2l3,2)42(2,3ll 

2C(l,?)  •  -(2l7,l)A2(3,3)-2(3,l)*2<2,3n 

2C(  1,31  •  FI2I2,  l)92(3,?)-£(3tl)P2(2,2n 


2CI2,ll  ■  -I2I1,?)PZI3,3)-2(3,2)42(1,3)) 
7CI2,2)  «  4(21 1,1) 42(3,31*213,1197(1,31) 
7Cf7,3)  -  -(211, 1)*2(3, 21*2(3, 1)42(1, 21) 


2C<3,!)  •  4I2M, 2)42(2, 31-2(2, 2)92f  1,31) 

ZCI1,2)  •  -(2(  I, 1)4212,3)^2(2, 1192(1, 311 


ACP04960 
ArP04‘>70 
ACP049R0 
ArR04990 
ACROSOOO 
ATROSOIO 
4ra05020 
Ar3  0‘’O30 
ACR0S')40 
ACPOSOAO 
4CR 95040 
ACROSOrO 
ACR050AC 
Ar.«05090 
ACROSl  T»0 
ACROSt 10 
ACROSI 70 
ACR0S130 
ACROS140 
ACR05150 
ACR0S140 
ATPOAlTf) 
ACRnStRO 
ArRC5190 
ACR05700 
AC»05210 
AfR0S?2O 
ArR052)0 
A(;R  05240 
ACR05250 
ACR05250 
Aril052T0 
ACROSZPO 
AC® 05290 
ACB 05300 
AC»05310 
ACR05370 
ACR05330 

ACP0534r> 

ACR05350 
ACR05360 
ACR053T0 
ACR05350 
ACR05390 
ACR05400 
ACR 05410 
ACR 05420 
ACR 05430 
ACRQ5440 
ACR 05450 
ACR05440 
ACR05470 
ACR0S4R0 
ACR05490 
ACR05500 


n  n 
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Fite:  *c*i  Fo»m*N  » 


OHIO  UNIVFItSITV  ElECT)(0>l*GNEriC  COHP«T*fl  IL  I T  lES  L«B( 


EC(3>3I  •  *l2(t.tl*Z<2.?i-/l2.tl*;il.2>> 

4rR055lO 

c 

4CR0SS20 

c 

TO  OBTAIN  ADJOINT  NATRIX  IZAI 

4C4 05^30 

00  I  J  •  1.3 

Af:R0^^40 

00  1  1  >  1.3 

ACRO55S0 

1 

ZAII.JI  >  ICIJ.II 

ArR05560 

c 

ArRf>9570 

c 

TO  OETERNINE  TMf  VALUE  OF  OETERNINANT  OF  IZI 

ACROSSRO 

OET  -  IZIl. 11*1712. ?I*;|3. 31-/13, 7IR/I2. 3)11 

ACROS^^O 

2  -1211.21*17  12,  11*71  3,3)-2(3,i)*2«2,3)ll 

Ar.»0540a 

3  *121 1.31*1 712. 1 1*713,21-713.11*2(7.21)1 

4CR0S4I  0 

c 

ACR0S620 

c 

TO  OBTAIN  INVERSE  MATRIX  OF  121) 

ACR0S630 

00  2  J  -  1,3 

Ar.»0S440 

on  2  1  >  1.3 

ACR09640 

2 

2111. J)  >  2A(I.J)/0ET 

ACR0S&60 

RETURN 

ACR0S670 

ENO 

ACR0^6R0 

c 

ACROS640 

SUBROUTINE  CONST  IV, F, ALPHA) 

ACR 05700 

c 

ALPHAS  OEPENO  ON  FREOUENCV  IN  NH7  ANO  GROUND  RESISTIVITY  IN 

ACR05710 

c 

OHN.M.  UNOER  HEAVY  RAIN  CONOITION  GROUND  RESISTIVITY  IS  TS.O 

O.H.  ACR05770 

niNENSION  ALPHA  (lOl.VLINEI lOl.Bf TAIJOI 

ACR05730 

F  •  F/i.aaE3 

ACROS740 

c 

TO  find  the  VOLTAGE  CLASS 

ACR057S0 

VLINEIll  -  362.00 

ACR05760 

VLINEI2I  •  SSO.OO 

ACR0577f> 

VIINE(3)  •  BOO. 00 

ACR057R0 

VIINEIAI  •  1200.0 

ACR05740 

00  I  1  •  1,6 

ACROSSOO 

OELTA  •  VLlNEin  ~  V 

ACR05A10 

IFIOFLTA.GE.IO.OO)  GO  TO  2 

ACR 05820 

1 

continue 

ACR'^^AAO 

2 

IF( I.EO.l 1  GO  TO  3 

ACR05R40 

IFI I.F0.2I  GO  TO  6 

Ar.R05850 

IF( I.Ea,3l  CO  TO  S 

ACR05960 

IFII.EQ.A)  GO  TO  6 

4CR05870 

) 

BETAIl)  •  8. OOF-6 

At* 05880 

BETA(2)  •  60. OOF-6 

ACR058<>0 

"ETAI3)  =  3S0.00E-6 

ACR05900 

CO  TO  7 

ACR05910 

4 

BFrAlll  -  9.10F-6 

ACR059?0 

BET1I2I  »  TO.OOE-6 

ArR059)n 

BETAI3I  «  3S0.00E-6 

ArR05940 

GO  TO  7 

ACR05950 

5 

BETA! 11  •  tO.OOE-A 

ACR 05980 

BETA! 2)  ■  70.00E-6 

ACR05970 

BETAI3)  -  3S0.00E-6 

ACR 05980 

CO  TO  7 

ACR 05990 

6 

BETAIl)  •  I0.60E-6 

ACR 06000 

BFTAI2)  «  86.00E'6 

ACR 06010 

BFTAI3I  ■  3S0.00E-6 

ACR06020 

GO  TO  7 

ACR 06050 

7 

CONTINUE 

ArR06040 

00  8  1  >  1,3 

ACR06050 

ac»i 


FORTRAN  A 


OHIO  UNtVbRSirv  fLECTROHAGNFTir  r.n*<PA  TAB  1 1  !  T  !  f  ^  L  ABO 


ALPHAM)  *  |f**('.R|ASORTIO*  ?SOI»R€TA<  II 
WRTTF 

FORHATI'O' ,nx, 'VALUFS  OF  ALPHAS*! 

WRrrF  (6»10»  (  I . AL'HAf I  1 ti^I •)! 

FOrtMATUOK,*  ALPHA  (SH#*!  -••GtO.1l 

RFTURN 

FNO 

SUBROUTINE  IPHASf  i  •< .  K  ♦  KK  ,kkk  •  IC?  I  •  IT  ^1  •  I C??,  IC5? .  IC2 1,  1C  51  •  1 P ! 
RFAL  30,10»,IC?U?0,?OI*IC5U?0.?0)tlPI20»?OI 

REAL  IC22C 20,20  »* 1C S?<?0,20) 

REAL  IC?1{20.20!,ir51(20,20l 
REAL  HAtlO.IO) 

on  I  J  =  Itl 
no  1  I  *  1,1 
MA( I , ji  s  Mf j, n 
no  2  I  *  1,1 

iFfK.Eo.n  ic(i,n  =  ic2uun 
if(k.eo.2)  tc<i,i)  «  trslif.ti 

IFIKX.fO.t)  *  IC22(f,n 

|F(<K,E0.2)  =  irS?Mt!! 

If{<KK.FO,ll  TC(t,t)  - 
|F(KKK.F0.?)  -  trsif!,!! 

CONTINUE 

no  1  I  *  1 . 1 

00  1  J  *  U1 

|P( f • J!  -  NA( t, J I «f C  t  I ,  t  I 

OF  TURN 

END 

SUPROUTINE  NOISE  ( K , XLOC ,YtOC tH*0»I P2 U I P? !. I P22, I P52f I p2 1, I P51, 
EE  I 

REAL  |0l< 20.20) , I P2T 20. 20 ) , I  PI 1 20.20! 
real  10211 30,301, 1051(30,101 

PEAL  (022120,201, !0S2( 20, 201.1 021(20,20 1.1053120, 20! 

OMENS  I  ON  EE  (30.  JO)  .  El(  20,20  I  ,E2(  20,20),E3f  20,20! 
on  1  J  «  1,3 

on  I  I  *  1.1 

rF(<.FQ.l  )  GU  TO  10 
IEM.E0.2!  GO  TO  !  I 
MKI.J!  »  !P?M!,.n 
I02( I, JI  *  IP22( (,JI 
fPIII.JI  =  IP21f  f,  JI 
GO  TO  I 

!PM  I.JI  a  IP5M  l.JI 
r02( If Jt  *  IP52(  f ,  J) 

|01( l.JI  -  (051(  I,  J) 

continue 

X  =  XIOC 
Y  =  X*0 
7  ^  t-O 
w  =  Yinr.*M 

U  »  YLOr-H 


ACROftO^O 
AFRO#, 070 
ACR‘y#»0R0 
ArRO^ORO 
ACR06I 00 
AC»fj6l  10 
AroO5l20 
ACROBlIO 
ArnOM^O 
ACRO^ISO 
ACR 06160 
ArR06  1  70 
ACRn6! «0 
ACR061R0 
AC R 06 200 
ACR06  2  1  0 
ArR06?20 
Ar«0*>210 
ArP06240 
Af  006250 
ArR06260 
Aro06270 
irR‘j6280 
ACR06200 
Arft06100 
Ar.ft063lO 
ACR  06120 
ACR06110 
ACR061A0 
Ar»06150 
Aro06360 
ACROMTO 
Arao6i50 
Af:R06lR0 
ACR064nO 
Kf  004410 
ACR06420 
ACR 064 10 
Ar»06440 
ArR06450 
ACR06460 
ACR064T0 
ACR 06450 
ACR06490 
ACR 06500 
ACR06SV0 
ACR06520 
ArR06530 
ACR06S40 
ACR06550 
ACR 06560 
ArR065T0 
ACRO6580 
ACR 06590 
ACR06600 
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At«l  FQHr«A»  OHIO  UNIVEH^ITY  (L  f  C  TftOflACNF  T IC  COHPATA'H  R  I  T  IFS  tARf 


FYI  »  A0*00*t  (W/l  Y*Y  ♦  W*WU  -  |U/IY*V  ♦  U^Ufll 

FXl  s  60*0n*Y*<  (  1«00/1Y*Y  ♦  U*(l||  >  <UOO/tV#Y  »  W»W)II 

FF(ttt)  «  SQPT(FY1«FY1  ♦  FX|«FXn 

FV?  *  &0,00*f  I  W/(  X*X  ♦  W*UII  -  IU/IX*X  * 

Fx?  «  6o.ooFx«n  1  .io/ix*x  «  u«uti  -  ii.oo/(x*x  *  wAwm 
FF<?#ll  ■  ^0XTIFY2*FY2  ♦  FX?FFX.t 

FVl  ■  60.00*1  (W/I2*Z  ♦  U«WII  -  lii/<2*2  * 

Ftf3  «  60.00*2*((1.00/f2*2  ♦  U*Ult  -  11.00/12*2  *  W*Htl) 
FFOtll  *  S0PT<FY5*FV1  ♦  FX)*FX1I 

TO  multiply  TmF  MOOAt  phasf  cuppfnts  mtth  The  field  factors 

FF!I*J>  TO  ORTAI^  THE  PI  NOISE  ifVEL  IN  NfCPOY/N 
00  2  1  «  !•) 

J  s  \ 

Ell  I . Jt  »  0.0 
on  ?  INDEX  *  1,3 

EKT.JI  *  EIM.JI  ♦  IPIIf,  rNOEXiPFFfiNOEV,  J) 

DO  T  I  s  1,3 
J  a  I 

E?fltJ)  •  0.0 
00  3  INDEX  »  1 ,3 

E?<T«J>  «  E7I1.JI  ♦  lP?n,INOEX|*FFIINOEX,Jl 

00  4  I  •  1,3 
J  *  I 

E3I  l,JI  -  0.0 
DO  4  INDEX  •  1,3 

F3nfJ1  «  E3U«J»  ♦  tP3U  •INOEXiPFff  index, J| 

E6I  a  0.0 
J  a  1 

DO  5  I  »  1,3 

^El  •  FFl  ♦  SORTIEll  I.JIAFK  I.JM 

'^E?  ^  0-0 

.1  =  I 

00  6  I  «  1,3 

FF?  =  Ff?  ♦  S0PT(F?ri,JI*E2f l«J»» 

FF^  ^  0-0 
J  *  I 

DO  T  I  *  1,3 

FFl  ■  EF3  ♦  SORT(E3(I,JI*E3f  UJtl 

! TOT  *  SOPT(EEl*EEI  ♦  EE2«FE2  ♦  EE3*EE3I 
FF  a  ?0.00*ALOGI0IETOT» 
i  ptupn 

FNO 

«;U4»0UT|NE  SIGNAL  (INC  ,K  ,  XP,Ht ,  H2,  Pf  ,  EPS  I  ,S  f  GMA  ,E1IP,C  ,  FFSI 
PEBl  LAN0OA 

COMPLEX  eta,oelta,o,rho,frf,epfc,t,tt,a 


ACA066I0 

ArP06620 

BCRO6630 

ACA0664Q 

acroaaso 

ACR06660 

acroaato 

BCX066A0 

ACP06600 

ACROATOO 

ftCA067IO 

ArR06  220 

ACX06T30 

ACRO^TAO 

4r.P067S0 

ACR06760 

BrRO*‘T70 

4CP06780 

BCft067P0 

ACA06R00 

AtROAAlO 

ACA06820 

ACRD6810 

ACR06840 

aC*068S0 

ACP06860 

ACX06870 

ACR06880 

ACP068N0 

ACR 06000 

aCXOAOlO 

ACR06O20 

ACRD6O30 

ACR06940 

ACB 06950 

ACRO6960 

ACR 06970 

ACRO6980 

AC*06«9D 

ACPOTOOO 

aCROTOlO 

ArRO7O20 

ACR0T030 

ACR070AO 

ACR07050 

ACR 07060 

ACR07070 

ACR 07080 

ACR0709Q 

ACR07100 

ACRO71I0 

aCROTlZO 

ACR07130 

ACR07140 

ACR07150 


FIL^:  ftCRl  FORTRAN  ^  OHIO  U?^IVf»SITV  f  L  FC  TR  OHAGSF  T  I C  C  OH®  A  T  AB  I  L  I  T  |  E  S  I  ABO 


C 


C 


C 


Cn^-PIEX 

1F(K.FQ«11  F  =  ?t)0.00 
IFtK.EQ.?)  F  *  SOO.OO 

EE  *  -UflOOE7*SIGHA/F 
ETA  «  CMPtXIEPSf.EFI 
delta  »  CSOftTIETA-1 .OOt/ETA 
LAMBDA  «  C/<F®l.OOFl) 

CK  »  2.004P1/LAMB0A 

«  *  XP*l«5?.00  ♦  FLOAT(INC|PO.T049 
00  »  SOBf(«*P  ♦ 

P«  *  SORT<P|*CK40n/2.00) 

00  *  PI /4. 00 
Q  «  CMPLXIO.OO.QOI 
PHO  s  rEX®(0)*RR 

r  *  cfxpro»ASORr(c»f®fiO/?.oo)FOELrAAii.oo«‘efHUH2>/fOEtTA*oo>n 

11  *  z 

12  *  0*1.0) 

ZT  *  /(5,0*<>.0*U0) 

74  a  |;•♦7^/(T,0*T.O*7•0♦U0) 

ERF  *  <Zl-Z?*Z^-Z4)*( ?.OOZSQ«TIP| U 
ERFC  *  1,000-ERF 

T  a 

TT  a  C€XP( T) 

A  *  i,oo-<RHn*neLTA*TTAeRfci 

ES  «  <0,4B7*$0RnERP|/(R)|PCAflS<AlPU00f6 

EF$  *  70,00*ALOGl0< E$l 

RETURN 

END 


ACROTIBO 

ACROTITO 

ACROTIRQ 

ACR07190 

ACR0T700 

ACR072I0 

ArR07??0 

ACR0T2T0 

ACR07240 

ArR0T250 

ACR07Z60 

ACR07270 

ACR07?A0 

ACR07290 

ACROTTOO 

AraoTilO 

Arponzo 

ATROTiXO 
ACftOTT40 
AFR  073S0 
ACR073A0 
ACRO73T0 
ACR073R0 
Ar.R0/3OO 
AfROTAOD 
ACROTAIO 
ACR074?0 
ACR07430 
ACR07A40 
ACR074BO 
ACR07440 
ACR07470 
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MM-  SIGMM  FORTRAN  A  flH  I  0  t>N|  VfRS  I  TV  f  L  fCTROHAGNe  tlC  Cfl^RAT  AR  1 1  I  T  f  F  S  LA«f  . 


r  TO  rUMOUTF  THF  SIGNAL  STRFNGTH  OF  AN  NOR  TRANSNITTFR 

r  GROUNH  CONDOCTIVf TV  IS  10  M^hO/M 
C  RELATIVF  PERMITTIVITY  IS  10 

r.  FFFFCIIVF  RAOIATEO  POMFP  IS  1  Kw 

C 

OIMFNSinN  YI200I  ,K(  ?00»  .Y1  (  121  ,XIU?) 

DATA  Yl /'S*  ,  •  !•  ,*r.«  ,»N*,  •  A*  , ‘L*  •• 

DATA  KL /'O'  iM  •  S*  f  •  T*  ,  • 

COMPl FX  FTAtOFl«Q,PQ,2,FRf , FRFC,T,TT*A 
COMPiFx  71, 22,73*24 
C 

hat  A  ‘M,F,Sir.,P,C/3.  1416,  10.0,0.010,  l0.0Q,3.0Ee/ 

r 

WRITE  (6,201 

20  FORMAT!  •  !•  ,  T20  ,  •  f)  1ST  I  NM  I  •  ,T  IS,  *FSOR<  200>«  ,TSO,  •  ESOBI  5001  •  > 

C 

C 

XF  a  100. 0 
XI  =  0,0 
NP  «  100,0 

OX  S  (  XF^Yf l/FLOATINP) 

0  s  X( 

00  I  I  »  1,N“ 

0  »  n*ox 
J  »  !  ♦  100 
HI  •  40.0*0.3049 
H  a  0000,0 
H?  *  H*0,3049 

c 

C  FOR  F PE  0  Of  200  KHZ 

c 

E  »  200.0 

FF  *  -l,BE  7*SIG/F 

eta  »  CMPLX(E,FF) 

c 

liFL  a  CSQRT<6TA-l.OI/FTA 

r, 

WL  =  C/(F*l.0F3l 
CK  a  ?.0*P(/Wt 

!>0  =  SOPTno*l«S2,0»**2  ♦  UU-H?I**2I 
f»»  s  S0RT(PtECK*OO/2,0l 
CO  *  PI /4.0 
0  ^  rMPLXIO. 0,001 
pn  «  CEXP(0)«PR 

c 

7  -  f  FXP(  01  •SOP  T(rK*nn/2.ni  POFL*!  l.0FTfHI^H2l  7f  OFlPOnn  I 

c 

C  TO  COMPUTE  the  FRFCIZI 

c 

7  1  »  2 

72  »  <7PP3)/(3, 0*1.0) 

73  *  I7**5I/(5.0*2.0*I,0) 

74  «  I2**7) /(7.0*3«0*?.0*l .0) 

r 

FRF  =  I  71  •72*21-24t*I 2.0/S0Rr(P| U 


SIGOOOlO 
SIG00020 
S  IGOOOTO 
Sir, 00040 
SIGOOOSO 
SIGOOOSO 
f>  IGOOO  70 
SIGOon«o 
S IG00090 
SI  0001 00 
S  IGOOl 10 
SIG00120 
SIG0OI3O 
S  IG00140 
SIGCOISO 
51000160 
SlGOOl 70 
SIGOOl «0 
S1G00190 
SIG00200 
SIG00210 
S IG00220 
SIG00230 
SIG00240 
SIG00250 
Sir.00260 

SIG0''270 
SI0002B0 
S  |G00?P(7 
SIC00300 
SIG00310 
SIG003?0 
$100-^330 
51000340 
SIS003S0 
SIG00A60 
5IG03370 
SIG003fl0 
SIG0039t> 
SIG00400 
SIG00410 
SIG00420 
SIG0043O 
SfG00440 
SIG00450 
SIG00460 
51000470 
SIG004R0 
$1000490 
SIG00500 
SIGOOSIO 
SIGOO570 
SIG00530 
SIG00540 
SIG00S50 


wwuw 


I 


1S2 

i 


PIICI  SI*«ML  V(Wr«*N  *  OHIO  UNIVEKSITV  EICC  TUnHAGNETfC  COHRATABIl  ITIFS  l*BO 


nvr  •  i.a-EBF 
e 

T  •  IBI 
TT  ■  fEXFin 

t  m  |.e-IIIO*OEl*Tr«E«FCI 

e 

II  «  (•.MT*SO«riFI/in*iaB?.OII*CAB$UI*UOE6 


10*  F*M  OF  BOO  KHZ 

Ft  •  F  •  100.0 
Fll  ■  -1.0fT*SI6/Ft 
FT*  ■  CHFLIII.CFIt 

Dll  ■  eSO«T<fT*-t.OI/ET« 

HI  •  e/(F|*I.OEll 

CR  «  F.OFFI/IH. 

00  ■  $0*THO*I*1Z.OI**2  *  (H|-HZ|«*2t 
*•  •  SO*TI*I*CK*00/Z.0l 
00  •  FIM.O 
0  ■  CMFlXfO.AfOOl 
*0  «  C««F(OI*** 

Z  ■  r.fV*IOt*10»l(CK*O9/2.0>FnEt*lt.O»IIHt»H2)/IDFt*Ot»M 
Zl  ■  Z 

It  •  IZ««1I/(1.0*1.0I 
Zl  •  IZFFlI/ll.OFZ.OFI.OI 
ZB  a  IZ**TI/(T.0«3.O*Z.O«I.O) 

F*F  •  IZl>tZ*Zl-ZBI*(?.0/$Q«T<Ft|| 
f*Fr  ■  t.O-EBF 

T  ■  ZFZ 
TT  a  CEXFITI 

•  a  |.0>fRO*OEL*TTaE*FCI 

FZ  a  l9.BHT*SO«T|FI/(Da|BlZ.OII«C«BSI*t*I.OE6 

FMI  a  ZO.OFEIOCIOIEII 
Kill  a  0 

VIJI  a  ZO.OFUOCIOIfZI 
KIJI  a  0 

HRITF  10,101  Kill, Till, T(JI 
to  Fa«H*T|lOX,lFll.ll 

I  FONTINHC 

rou  HFiOTix,v,  100, z,n, Til 

WITF  10,091  M 

99  Fm9*ri*0',t0X,*REC.  HEIGHT  a«,F|0.9> 

ITOF 

FNO 


$  IGOOBBO 
IIGOOBTO 
IIGOOBBO 

Sfcoosoo 

SIGOOBOO 
SIGOOBIO 
SIG006Z0 
S 1000630 
B IG00660 
SIG006S0 
SIG00660 
S 1000670 
StC006B0 
SIG00690 
S IG00700 
SIG007I0 
SIG00T70 
11600730 
SI600760 
SI6007SO 
$1600760 
$1600770 
$1600700 
$1600790 
$1600*00 
$1600910 
SIGOOBZO 
$ IGOOBIO 
$1600840 
$ IGOOBBO 
$1600860 
$1600670 
BIGOOBBO 
$1600*90 
$1600900 
$1600910 
$1600970 
$1600930 
$16009*0 
$1600930 
$1600960 
$1600970 
$1600960 
$1600990 
$1601000 
$1601010 
$1601020 
$1601030 
$1601040 
$1601030 
$1601060 
SICOIOTO 
SIGOIOBO 
$1601090 
31601100 


n  n  n  n  n  n 
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FlieJ  OCRI  FORTRAN  A  OHIO  UNIVFRStTY  f LFC TR OH AGNE T 1C  COHOATAB IL I T  IFS  IAB( 


C  TO  COMOUTF  the  »I  NOISE  OF  THE  OC-t f NE  WITH  THE  POINT 

C  OF  ORSERVANCE  at  any  tOCATION. 

c 

C  CONSIDER  ONtY  THE  POSITIVE  POLE  OF  THE  BI-POLAR  AS 

C  THE  NEGATIVE  POLE  GIVES  A  NFCIIG|PLE  P|  LEVEL. 

C 

niMENStON  YI200)  ,X(  ?OOI  ♦YU  t?)  •  VI 1 121 

DATA  YL/*  •f*R»,M%*  •**L«**E««»V***E*#*L*f* 

DATA  XL/*0*t*I*, ‘N*.*  *. "F* E T* / 

C 

DATA  oi,V,F/T, 1416, /,00. 0,200.0/ 

DATA  NSUA,«A0#H,H,S/4,U52S,45.r,4R.RS.'T6.T3/ 

C 

ENORl (Ell  *  El 
FN0P?(E2)  =  E2 
FNORIIETI  s  E3 
EN0R4(E4I  «  E4 
C 

C  TO  DPTERRfNF  ThF  MAXIMUM  SURFACE  GRADIENT. 

C 

R  R/(?.0*SfN(PI/FLOAT(NSUflm 
A  »  1.0/F10AT|NSI)R) 

N  «  NSUR-l 

P  »  (EL^ATfNSUA I «RA04(R 4«NI |P*A 
0  *  SORT(  ( I2,0AH/S|P42»  ♦  U0> 

OG  •  ?,0RM*O. 3040*100. 0/<P*0l 

G  •  (l.O  ♦  (ELOATINI*RAD/R||/(ELOATfNSURt*RAn*AtOG(GGU 

C 

GMAX  «  0*V 

C 

WPfTC  (6,11  F,V*NSUR»RAD*B«H«S»G,6NAX 
I  FORMAT  (M  •  •?0X,  •FREOTKHZ  la*  tE 1 0. 3//2 1 X»  •  VOL  T  AGE ( t-C-KV  I •*  » 

6  Fl0.3//21x**NSUPs*.l4//2lX**RA0lUS  OF  SUR(Cm«*t 

6  FI0.4//2I X, •«  (CM|a*,F|0.3//2|X**AVE  HEIGHT  fFT»s* 

6  .F10.3//21X, •POLE-POLE  (FTI  »• . FI 0. 3//2 1  X, 

6  »G  (KV/CM/XV)  *• .F10.4//21X, 'GNA*  (KV/CNI  a*, 

6  F10.4//1 


TO  rf'»<PUTF  THF  R|  AT  GROUND  LEVEL  AND  FROM  THERF  TO 
DETERMINE  TMf  TOTAL  CURRENT  IN  THE  CONOUCTOR..  THEN 
USE  THE  SU«R0UTINE  HEIGHT  TO  COMOUTE  THE  PI  AT  ANY  LOCATION. 


WRITE  (At4l 

4  FORMATCTl  W*  OTST(FT  I  •  »T?4  •  •  HP(  O.OETt  •  t  T36t  *HP  (  ?50FT  I  •  , 

6  r49,*HP( AOOETI •♦T60,*HPf tOOOFTI*/! 

C 


XF 

-  900. 

,0 

XI 

*  0.0 

NP 

*  90 

nx 

*  (XF- 

«n/fLn4HNP» 

0  > 

‘  XI 

00 

2  I  » 

l.NP 

D  *  0  ♦  OX 


OCROODIO 
OCR  00070 
DCR00030 
DCR00040 
DfROOOSO 
0CR00060 
OCROOOTO 
OrROOORO 
DCROnORO 
DCROOlOO 
OCROOllO 
0CR00I20 
OCR  00  I  30 
or  ROD I  40 
DCROOISO 

0r»00160 

OCR  001 70 

OCRDOIAO 

0CRO0I90 

0C»00?00 

DrR002l0 

0CR00270 

0CR00230 

OCR  00740 

nf.ROO?50 

0CR0O760 

DrRD07T0 

0CR002B0 

OCR  00240 

OCROOBOO 

nCR003i0 

DCR00370 

OCR 00 3 30 

Or.R00  340 

nrnooiSO 

0CR00360 

OCROOSTO 

DCR003R0 

OCR 00 390 

nCR00400 

OCR 004 10 

0CR00420 

OCR 004 3 9 

nrR00440 

DCR004SQ 

OCR  00460 

0CRO04T0 

OCRO04R0 

0r»00490 

OCROOSOO 

OCR  009 10 

OCR  009  70 

OCR 00930 

OCR  00940 

OCR00990 
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nCRI  FOKTRAN  A  OHIO  UNIVERSITV  ELEC TAOWACIET 1C 

CnNPAr«flILIT(FS  LAQf 

J  •  I  »  30 

DCROOS60 

K  •  J  *  30 

OCR00S7O 

L  •  A  ♦  30 

oc»oos«o 

c 

DrP00S90 

c 

DC ft 00600 

C>AT  a  CRAAFM.OO 

OCft 00610 

OAAO  a  3Q«TIH*H*0*nl*0.3043 

OCft006?0 

c 

Drft006^0 

E  •  ElA.OAALOGtOIGAAn  -  7TA.O*ALaCtOtCIIATI*AiaGIOICPAT  1 

OCft00640 

?  *  40.0*AL0G10(AA0I  *  2T.o*AinGtOI83A.0/F| 

0rft00660 

3  *  «0.0*ALaGIOM0.3/0«AOI 

OCft00660 

c 

Drft00670 

c 

TO  COR»ulf  THE  «!  MOJSE  AT  OIFFEAENT  HEIGHT  AMOVE  GROUND 

r>fft00660 

c 

AT  THE  HEIGHT  OF  0  FT. 

Or» 00690 

HF  a  0.0 

orftooToo 

CALL  HEIGHT  lE.H.HO.O.FI.EHl 

OCftOOTlO 

El  a  EH 

OCft007?0 

c 

OCRO071O 

r 

AT  THE  HEIGHT  OF  300  FT. 

DCft00740 

HR  a  230.0 

0rft00760 

CALL  HEIGHT  |E.H,HF,0,»l,EHI 

0Cft:>0760 

f7  -  EH 

f>fft00770 

c 

r>r»007ft0 

c 

AT  THE  HEIGHT  OF  1000  FT. 

OCftOOT90 

M»  a  300.0 

or.ftOO%oo 

CALL  HFIGHT  (  E.H.HF.O.FI .EHI 

OCftOOftlO 

G3  a  EH 

ocftooA;»o 

c 

nrftooeio 

c 

AT  THE  HEIGHT  OF  3000  FT. 

0rft00R40 

HF  a  1000.0 

OCftOOftSO 

CALL  HEIGHT  1  E>H,HF.n,F| , EH) 

orftoo«60 

E4  a  eh 

OCftOOftTO 

c 

OCftOOft«0 

c 

orftooft9o 

Till  a  ENDOIIEI 1 

DCft00<»00 

XIII  a  n 

r>r«oo9io 

FIJI  a  FNOO?IC?I 

or«oo97r) 

XIJI  a  0 

DCft 009^0 

TIKI  a  ENORTIEXI 

OCR  00940 

XIXI  a  0 

OCft 00950 

VILI  •  ENORAIEAI 

DCft 00960 

XILI  a  0 

t>Cft00970 

c 

OCft009ftO 

c 

or ft  00990 

NRITF  lAtSI  XIII. Till. TUI, TIKI. Till 

OCftOIOOO 

3 

F0RNAT|3X,3F17.3I 

OCftOlOlO 

? 

CONTINUE 

DCft 01 070 

CALL  HNFIOT  IX.T,30.A.Xl.TLI 

0C»0I050 

STOF 

DCftOI040 

END 

OC«01050 

c 

DCft 01 060 

r 

TO  FUT  IN  SURROUTINE  HEIGHT 

ncftoiOTo 

c 

OCftOtOftO 

SURROUTINE  HEIGHT  1 E.H.HF.O.FI.EHI 

or.ft  01090 

c 

OCftOllOO 

1S5 


FILE:  nC«I  FHRIRftN  A  OHIO  ON|VFf»^ITV  f  I FC  TWO-^ACNF  T  1C  COHP  AT4B  1 1 1  T  I  F  S  LAB^ 


c 

f  IS  IN  Oft/Hirv/H 

DCR0I110 

c 

FF  IS  IN  HICV/H 

0rR01l20 

c 

OCROI 130 

»AT  =  F/?0.0 

OCROIIAO 

FF  s  10.0**BAT 

O^HOl ISO 

c 

OCR01160 

c 

TO  GFT  CUBBFNT. 

OCROI 170 

ruftB  »  FF*nH»0.304<JM»2  ♦  1  0*0.  3049l*#2) /I  60. 0*2. 0* 

OCROI IHO 

6  M*0.306<?» 

0CR0U90 

c 

nCRO|?n'i 

c 

NOW  USE  CURRENT  TO  COMPUTE  H-F|flO  AT  ANY  OOINT. 

or»ot  71 0 

c 

r)CR0l?70 

p  s  *0. 3049 

0r.O0l230 

0  s  (HR-H»*0.3049 

ncRoi?40 

on  s  0»0.1049 

nc»oi2so 

c 

OCROI 260 

HT  a  (CllRR/<2.0*PfM*(<P/f0f)^0I>  ♦  PRPJI  -  10/100*00  *>  0*0n) 

ocRni?7o 

HY  a  (rURR*00/l?.0*PIM*UUO/<00*00  ♦  0*0))  -  (1.0/ 

nrftoi 

6  (00*00  ♦  o«PI ) I 

nrRoi?Ro 

c 

OCROI 300 

FK  *  taO.OPOlPHX 

nrRoi3!o 

6Y  *  120.0*PI*HV 

DCR01320 

FHP  a  SORTIFX*?*  ♦  «y*eY) 

nCR0l330 

FH  *  20.0*ALOGIO( «HO) 

OCROI 340 

RETURN 

OCROI 3S0 

END 

O'-ROI  360 

156 


nif:  PJIWICH  FOWTOftN  K  OHIO  UNiVfRSITV  f  LFC  TR  0*<  ACNF  T  IC  COHRA  T4 1 1  1  F  IF  ^  LAB^ 


I^'PlICfT  RF4l*PIA-H,0-;,$) 

Rf  Al «4  XM, XAR ( tOO 1 , YARI  1001 tXARL<  t2l •V4RLC 12) 
f EP?,FP3,FTA,CAH,vll,/lt*7S 
rnM0|^x«i6  Fux.Fuv.fu? 

f  ;:MPtcx*l6  rpo^.FPTSff  TP^,FTTS*f  X*F  Y,€  r  ♦U<*LO 

tOMPt  rx*l6  r  (400S  »  ,C  »FP|<»0»  •Fri90)  .FP«>|90)  TTI90) 

PI'AfNSinN  111901 , 12  1901  ,  f  31901  ♦  |A<90) 

rCl“OlFX*l6  com  991  ,SGO|  99)«rG(  19B)«VCI19B1  tZLOl  1991 

OIMFUSION  0(991,1 A|99),|9|99|,|sr|99),^0l99,4),N0l99) 

OIHENSION  XIIIOI.YI  1101  101 

DATA  PI ,TP/3« I  4159265 398979,6,?ft3185307|T9S9/ 

DATA  E0,U0/9.954E-12«1.2566F-6/ 

C 

c  FOR  THF  PLOT 

OATA  YA9t/«0*.M*,*S*,*T«,*  •,«  •N«.*  •,*F»,*T*/ 

DATA  YARL/*H*  ,«A«  ,*G*  *•*  «,M*#*N*.*  ♦,*  •.•n*,*R*,*  •/ 

C 

2  FOPMAT(1X«8F10«41 

3  FPPMA  T( |X,4F10«4/) 

4  FfJPMA  T  M  K,  1 15,5^10,4  1 

5  FOPRATIIHOI 

6  FPPMATC |X,6F10.4/ 1 

7  FrPIAFlAFlO.SI 

8  forma T| IX, l|4, 1315) 

9  fnRMAT|3X,««AX  «  • ,  1 5, 3X , • « I N  s  *t(5,1X,«N  «  •♦151 

10  FnRMAmx,2l5,2FlO.?l 

11  <^nRMATnRl5,51 
ICJ«90 
INH*99 

on  14  J«l, INH 
15  ISC(JI»0 

RFAOI 5,7}RM,FP2,SIG2»r02 
WRITE (6, 21 RM,ER2# $102,702 
®FAn(4,7|AM,CR»',€R3,SlG3*T03 
HR|TP|4, 2UM,Cpm,FR3iSIC3,  703 

RFA0(5,9) f RISC, IG ATM, INEARff SCAT, lWR,N6eN,NW,NP 
WP! TR 16,51 IR ISC, IGA  IN, INE AR , I  SC  AT , | WR tNGFNfNW, MP 
RFAOI 5,7)FMC,RHA,THA,PH| ,TM| ,pmS,THS 
WR«TF(6,2IFmc,PM4,ThA,PHI,THI,»hS,TKS 

on  2?  j»i,nm 

RFAOrS.^MAf  J),  |9(  Jl 
22  WP|TFI6,9IJ,IA(J),|QU) 
on  40  IsttMP 
RFAoi5,tiixcii,Y<  n,zui 
40  wPITf  (6,4M,xn),Y<  11,2(11 
BfA0(5,lllXP,YP,2P 
XXPsXP 
VVP»VP 
;2R*4. 

FH7*FM^t|.E6 

nMEGA^IPPFHZ 

IF(  S!  G2.tT*0.  IEP2-ER2PE040rHPlX(  U0000,-T07| 
|F|TO2.Lf,0.)6P2*0CMPlX(FR»4E0,-SIC2ZONEG4) 

IF(SIG3.tT.O. IFP3»ER3P€OPOCMPlX(l,OnOO,-T031 
IFf  Tn3,tT.O.»EP3*Or'iPtXlER3*EO,-'SlG3/OMEC4l 


RJL 00010 
RJL00020 
P JL0O03O 
R JL00040 
P  JtOOOSO 
R JL00060 

p jtoooro 
RJt 00090 
9 Jt  00090 
RJLOOlOO 
R JLOOl 10 
P JL00120 
R JLOOl 30 
R  Jt 00140 
R JL00150 
R JLOOIAO 
R  J100170 
R JLOOIRO 
R  JtOOI90 
P JL00200 
RJLO0210 
R  J100?>0 
P JtOO’30 
R  JLOo>40 
RJL 00240 
RJL 00260 
R JL002ro 
RJl 002A0 
RJL 00290 
PJL 00300 
RJLO03I0 
» JL0O320 
RJL00330 
RJL 00340 
RJl 00340 
R  Jl 00360 
R  JL  00370 
R  JL003R0 
P  JL00390 
RJl 00400 
R  J100410 
R  JL00420 
PJl 00430 
R J100440 
RJl 00450 
R J100460 
RJ100470 
R JL00450 
RJl 00490 
R  J100400 
R  JL005IO 
R JL00570 
RJ100530 
RJl 00540 
R J100550 


n 


]57 


f^ILF'.  «JlRICh  FDPTftiN  4  OH  10  »JNI  VfcRS  I  f  t  f  LFC  TR0H4GNE  T  I C  COMPA  Ti^l  L  I  T  I  r  s  LARO 


FT4*f.nS0RT(U0/FO^  » 

GAHanHEGAPCnSOft  » ( -UO*FP^» 

CALL  PLSORTf  I  A,  f  R  ,  !  1,  !?  .  H  •  JA  ,  JR  ,  MO,  NO,  KJS,NP,  N ,  ma  x  ,M  r  N  ,  fC  J*  r  NM| 

WPf TF(6*S» 

W0|TF«ft,9)H4X,M!N,N 
WPfTF (^,51 

!F<RAX.GT.4  .OR.  MIN.LT.l  .OR.  N.GT.trjJGO  TO  ROO 

1NT*A 

TI2*l 

no  60  J*l,NM 
VG(J»=<,0,.0> 

2Ln(ji=<.Of.oi 
J J*J>NH 

60  7L0(  JJI«UO,.Ot 

fF(NGFN.GT.01  VGI  NGFNT-f  I  - 

r  ALl  SGANTMA  ,!R,  INM,  INT,ISr.,  1 1  *12,  n*  J4,  JR,NO,  N,  NO.NM.NP 

2,AM,RM,C.CGO*fHH,0,FP2,6P3*eTA,FHZ,GAM,SG0,)f»  Y,Z#rL0,7S» 
TF(N,LF,0|Gn  TO  ROO 
fF(NGfN,L6.0IGn  TO  400 

CAIL  GANTin  A,I8,INM,  CWft*  I  If  I2»f  3#f  I2t  JA»  JR*HO,N,  NO*NH,  AM 

2tC.CJ,rGfCM«,n,6FF,r,AM,GGtrGOfSGOfVG*YUf?n*7inf7$) 

GO  TP  R67 
R63  no  ion  |«l,N 

1/2 

00  100  J*I fN 

I JaTO*J 

too  WR(TF(6«  tOI t f J«r< 1JI 
867  rONTTNUF 
WPlT£(6fS1 

H«|TFt6f3IFPF,GGfni 
200  TFUNFAR.IF.OIGO  TO  300 

I  OOP  TO  CAtrULATF  NFAR  ZONF  PATTERN 
INPUT  niMFNStONS  IN  FF6T 
no  \<i  FLeVATION  AROVE  0«nUNO 
tXO  IS  CLOSEST  DISTANCE  FROM  TRANSMITTER 
YYP  IS  fartmfst  distancf 


36S 

C 

C 

C 

c 


c 

c 

OR 


rPTSPl =10 
XMf NP(  sXXP 
xMfixei sVYP 

WOITF  (R,34S»  NPTSRLf  XMfNRLf  XMAXRL 
FORMAT  MSfTFlS.SI 

FOR  VAOIDUS  POINTS  OF  X,V,Z, 

TO  COMPUTE  values  AT  0|FFE»FNT  AtTITfOFS 
no  9QS  I7HT  »  If1 
7/P  »  ROO.O 

IFII7HT.FO.2l  77P  a  1000.00 
IF(I7MT.F0.3I  77P  *  ISOO.OO 


GO  TO  94 


R  JL 00*60 
R JLOOSTO 
R  JL00‘*3n 
R  Jl  00690 
R  JL00600 
ft JL 00610 
ft Jl 00620 
R  JL006TO 
P  Jl  00640 
R  M_006S0 
P  JLO'^-jSO 
R JL  00670 
R  Jl  006R0 
ft  Jl  00690 
RJl  00700 
ojl  ‘07\O 
ft J100720 
RJl  007TO 
P  jL0074n 
» JL  no7‘  0 

R JL007fO 
p JL00770 
ftJLOOTSO 
ft  JL00T90 
PJLOOftOO 
ftJLOOflIf) 
RJl 00R20 
R JLO0R10 
ft JLOOCAO 
R  JIOORAO 
RJL00R60 
RJL 00*^0 
RJl  00«P0 
R  Jl00ft90 
ft J100900 
R JL00910 
RJL  00920 
R  JL0n930 
ft  Jl 00940 
R JL009SO 
P  JL  00960 
RJl  00970 
RJL009R0 
R J100990 
P  JLOI  000 
PJlOIOlO 
R  JL01020 

R JLOIOAO 

ftJLOl 040 
ft JLOIOSO 
•  JL01060 
ft JtOlOTO 
ftJlOlOftO 
ft JL0I090 
ftJLOllOO 


CONTINUE 


FtL€t 

WI  i  FDftTRAN  A  HHin  UNIVfRSITY  €t€C TRHH ACNFT 1 C 

CO'^PATABILITIES  lAftf 

C 

PJLOinO 

C 

f-  -niNTS  ALONG  A  SFCTOR  OF  CONSTANT  RAQIUS* 

QJL0II20 

C 

P JLO! !T0 

ri'i  .*.7  1RJL«1  fNPTSRl 

ftJLOl 140 

P.‘  *  XitP4FtnATn«JL-ll*<VVP-XXP|/FtnAT<NPTSftL) 

PJLOllSO 

0''  ‘**A  '  fCUT»l  ,fiO 

PJL01160 

;p-  '7o/^.?«i 

BJLOUTO 

K  r  a  4  0!liSPnCOS{FLO4T(  ICUT-l  1  •PI  / 1  RO.O  1 

9  iLOl  ISO 

YP  r.  AO!  »JS*OS  INI  float*  ICUT- M  «P|  / 1  40*01 

RJL01190 

cn  19 

PJL0I200 

C 

PJL01210 

9T 

r  ON  M  N'l? 

•JL0I220 

C 

FOh  points  ALONG  THF  X-AXIS  ONLV 

PJL01230 

»l  -  A076.<.12*1*00  -  ?00.0 

RJL01240 

tr  6076.4l2»l*00  ♦  200.0 

ftJL0l?50 

M.»  50 

R JL01260 

SI 

«JtOl2TO 

r  «  f  t-r  »  1  1  /  •  flA  T(  '1 

ft JL0t240 

no  .  *( 

ftJl0l290 

on  mi  *  =  1 , nnp 

ftJLOllOO 

XP  »  <00  ♦  fioat<  1*  )<l/s.2Ai 

ftJLOniO 

YP  »  ).0 

ftJlOl3?0 

>p  ^ 

RJL01330 

00  T">  oq 

ft JL01340 

c 

ftJlul  450 

CONTIN'JF 

ftJL0n60 

C 

RJl.OnTO 

f 

FOR  PMNTS  tHF  Y-AXfS  ONLY 

•  JL0I3S0 

e 

RJL0!390 

YI  -  -‘-00.0 

ftJLOKOO 

VF  »  ‘1  >  ).0 

ftJlOlAlO 

NP  -s  |O0 

RJL01420 

NNP  a  im 

ftJlOlAIO 

OY  »  (  VC  -V  i  )  /f  1  OA  Tf  ^1>| 

ftJL01440 

00  «  YI 

ft J10I450 

on  PT  (  *  UNNP 

ftJL0l44O 

XP  s  »^^00.0/1,7«l 

ftJLOIATO 

VP  in.i  f  ri  OAfl  I-l  J*0Y)/3.2lll 

ftJLOlAftO 

rt'  '  .1 .  <  I  • 

RJ101490 

r,n 

ft JL01500 

c 

ftJlOlSlO 

96 

rnf.T  ■;!,* 

ftJLOlSEO 

c 

ftJtOISIO 

c 

pnpj.  ,  Mi  l.  (.  irt  ■  1  040  WITH  RFSPFCT  TO  TMF  OKC 

ftJLOISAO 

c 

T  P  A  N  •  -M  r  t  M- 

RJIOISSO 

c 

ftJLOlSOO 

X  r  »  s  .* !  , 

ftJlOISTO 

YF  s  i*.  M  7.  ! 

ftJLOISftO 

NP  ^  •  J 

ftJtOlS90 

NNP  Mi 

ftJtOUOO 

f>X  .  ,  1  »  /f  1  fM  T  f  ')P» 

ft JLOI6tO 

or  1 

ftJl.0l6?0 

0"  19 

ftJt0l630 

AN  '  . . .  (1  1  |♦^)X 

flJtOI640 

GP  A ‘ .  .  / . . 

ftJtOUSO 

f-ll  t:  FORTRAN  A 


OHIO  UNfVFftSIlV  fLfCTRO*<AGNFriC  COMRA  TAB  II I  T  I F  <;  LABO 


p..  =  -13324.09 

R JL0t660 

YM  s  GRAORXM  »  MR 

R.IL0I670 

KP  a  X'««/3.?Rl 

R JL016B0 

VP  =  VH/3.281 

RJL0t690 

;p  =  5.0/3.281 

R  JLOl  700 

r.n  TP  99 

RJLOl  710 

c 

R JLOl T>0 

<»1 

rONTINijE 

RJLOl  730 

c 

R  JLOl 740 

c 

“PINTS  ALONG  THF  AtCfSS  ROAO  WITH  RESPFfT  TO  THE  f.HH 

RJLOl 750 

c 

TRANSMITTER 

R JL0I760 

c 

RJL0I77O 

M  =  70715.0 

RJLOl 7R0 

*f  a  71215.0 

fl JL01790 

NP  a  ,>00 

RJLOIROO 

nnp  =  »ni 

RJLOl B!0 

OX  i  (XF-XD/HDAVINP) 

P JL01R20 

on  s  Ki 

RJl OJ  830 

no  84  I  a  l.NNP 

RJLOl  «40 

XM  a  nn  ♦  FLOATI I-l ) •ox 

R  J10I850 

G»AO  »  -0.655 

RJLOl 360 

BB  a  46475,83 

RJLOl 870 

VM  *  GRAORXM  ♦  BB 

RJLOl B80 

Xp  a  XM/3,28l 

R  JL01890 

yp  a  yh/3.281 

ft  JL0l9r>0 

7P  a  5.0/3,281 

R JLOi9!0 

GO  TO  99 

RJL01920 

c 

R JL01930 

«9 

rONTI  n'iF 

RJL01940 

c 

R  JL01950 

II 


92 

C 

c 

c 

c 

c 

c 

c 


(.MC  ONFlOt  I  A,  IB,  INM,  11,1?,  n,MO,N#NO,N'<,  A'^,CGO,SC'>,FTA,GA« 

J,0,t,Y,Z,XP,VP,2«»,ex,FV,F2,FUX,EUV*6U2l 
s  K0*9,?R1 
=  YO*3,?fll 
/■>  « 

WMfTF|s,3l)XO,VP,ZP 
^^RNAT  (3X,3tlO«!) 
v.Pr  TFI^,^6iFX,FV,  e?  ,etJX,FUVt€UT 
6^  ►  OP»^AT  (  31  ?X,20l?,4l  ) 

pATim  s  roABsiFxi/coABSiezi 

PATfOZ  *  COARSCEUXl/fOABSIFUZ) 

RATION  »  COABSIFUYI/rOABSieUZI 
WRITF  PATIO! .RATIO?, BATI09 

f UPNATI 19X,3F10«4) 


FI€LO  STOFNGTH  AOJUSTFO  FOR  IINI  9FACON  BY  »iUtTIPLYfNC  THE  VALUF 
OF  F9P  OSLO  IN  THF  PROOBAH  BY  A  FUOGF  FAfTOP.*,, 


FHR  uA'f  BEACHN  (  E«p*0.7M  I  THE  FIFIO  STRENGTH  AT  I  STATUTE 
F'!LF  IS  4933.3761  HICROV/H  OR  73.9629  OB/l  H|CV/N. 


r1>OGF^^.<.l0.327^ 

Ilf  LO*  ir  7»EUZI*lOflOOOO.O 
f  /»F7  •rjooooo. 
Euz«fu;*ioooooo. 


R  JV  0I9A0 
RJL0!9^0 
R  JL019RO 
R  J10IP90 
OJV  0?000 
P J102010 
R  JL02020 
R J102030 
RJL02O40 
R Jl 02050 
R J102060 
RJl 02070 
R JL020BO 
RJ102090 
OJL07IO0 
RJL021I0 
RJ102120 
R JL02I 30 
RJL02140 
RJL02I50 
RJt02160 
RJ102ITO 
RJL02180 
RJt 02190 
RJL 02200 


1  f  II  r  HP  1  3  AN 


fMii>  uNfvFp^irv  ri  Ff  T  1C  cnMPArft^Il  it  les  iabo 


’ ^  50 ,  *01 0 ~  1  nt  F  0  I'.i 

X  nAM',.;  t/  1 1 

« JLO??IO 

■^5n.*f)Ln'.ir)(ru  -.i 

*{  OAO«i(F‘jyi  \ 

R JL07?20 

•  .i5f)A-f  n  ’no 

0  JL02230 

-?n.*run  .iniF'f  ■ , 

X  nAft';(  F  i  1 

«  Jl 0?740 

,  '  ( y7i  li’ 

■  ,t!onn ,  'f',  <p,  7p,»iPt  n*T  ri  o« 

P  JL  0'’?50 

.  ,M  M  F 1  r. .  ^ 1  r  1 0 .  1 , 

10.  ;>,  1 F 10. 

ft  J10P540 
R  J102270 

l 

P  JLO??ftO 

.  • '  fl 

JLO7  5O0 
ft JL  02300 

'■  I  I  N'JF 

ft JL02310 

■  5  '  •.0',) 

ftJL  02370 

'  •  i'  .  NJf 

ojt 02330 

!'!  ■^no 

RJt02340 

•  .  !  N'jr 

P  JL023S’> 

.  '•(  QO*? 

P JL02360 

PJL07370 

1  1 .  •  I 

ft  Jl  0>3ft0 

;  iN'if 

P JL023O0 

■  f )  'j  .1  ■» 

ftJl07<.D0 

■  1  ■'  i  N<J' 

PJL  02410 

1  ( 

ft Jl 0?4?0 

■  •  1  r.A  INI  'J  ! 

P  Jl  02430 

•  ll 

P JL02440 

PJIOPAS'A 

'  -4A 

»JL02460 

•  —  r,  ►  ' .  '  [  ' . 

'  ,  ,  II  I3»  U?,Mn,N,  NO,N“.  AM 

ftJL0?47rt 

*>  T ,  r  .  -  • . ,  •  , 

CSP.rr.$Ttf*»*PT,F»P,€TT,FPPS*For$ 

ftJtO’AftO 

•  '  r  I  s  , 

•sr.  s  ;fsp,sc$t,^pph*sptm,stpm,sttm,  th 

ft  Jt074<i0 

^ .  / .  / , ;  i  n .  /  s ,  -  ’  4  ,  i  * 

ftJlO7S00 

i  rrji,,  ^  .T  , 

ft  Jin7si0 
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Appendix  C  Acronyms 
AC  Alternating  Current 

CMH  Designates  particular  NDB  near  Columbus,  Ohio 

dB  decibels 

DC  Direct  Current 

DKG  Designates  particular  NDB  near  Columbus,  Ohio 

ERP  Effective  Radiated  Power 

f  frequency.  Hertz 

KHz  Kilohertz 

kV  Ki lovolts 

MHz  megahertz 

NDB  Non-Directional  Beacon 

NM  nautical  miles 

RI  Radio  Interference 

s/n  signal  to  noise 

vs  versus 

relative  permittivity 
a  conductivity 


